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Abstract 

In this II, a probability to detect the neutrino produced in a high- 
energy pion decay is shown to receive the large finite-size correction. 
The neutrino interacts extremely weakly with matters and is described 
with a many-body wave function together with the pion and charged 
lepton. This wave function slowly approaches to an asymptotic form, 
which is probed by the neutrino. The whole process is described by 
an S-matrix of a finite-time interval, which couples with states of non- 
conserving kinetic energy, and the final states of a broad spectrum 
specific to a relativistic invariant system contribute to the positive 
semi-definite correction similar to diffraction of waves through a hole. 
This diffraction component for the neutrino becomes long range and 
stable under changes of the pion's energy. Moreover, it has a universal 
form that depends on the absolute neutrino mass. Thus a new method 
of measuring the absolute neutrino mass is suggested. 
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1 Introduction 



The method of computing the finite-size correction developed in I pQ is ap- 
plied to the probability to detect neutrinos in pion decay. Since a neutrino, 
charged lepton, and pion are described by a many-body wave function that 
follows Schrodinger equation, the kinetic energy at a finite time is not a con- 
stant and takes a wide range of values, as was shown in I. Consequently these 
waves include broad spectrum and show a diffraction phenomenon that is 
non-uniform in space-time. Since the speed of light is an accumulation point 
of the relativistic waves, a two-point correlation function has the light-cone 
singularity and the probability to detect the neutrino is subject to a finite- 
size correction. The neutrinos are very light and propagate with almost the 
light speed, hence the correction becomes unusual and depends, in fact, on 
absolute neutrino masses which are not found from oscillation experiments. 
Thus the correction becomes extremely large in its magnitude and size and 
has universal properties for the neutrino. 

Tritium beta decays j2] have been used for determining the absolute value 
but an existing upper bound for an effective electron neutrino mass-squared 
is of the order of 2 [eV 2 /c 4 ] and the mass is 0.3 — 1.3 [eV/c 2 ] from cosmological 
observations [3j. In these neutrino experiments, higher precision and more 
statistics have been achieved and will be improved more. 

Weak decay processes have been studied using an S-matrix of plane waves 
with asymptotic boundary conditions, where particles in the initial and final 
states are regarded as free waves without correlations jH El [6j [7J |8] . A decay 
rate, average life time, and various distributions of charged leptons have been 
computed, and perfect agreements with experiments have been obtained [9]. 
These facts proved that the standard theory of electro-weak interaction is 
correct. 

Neutrinos have almost the light speed and are detected at a distant po- 
sition, hence they are similar to an electrostatic potential of a moving body, 
which has a finite-size correction in a form of a retarded potential. Due to its 
extremely small mass, the probability to detect the neutrino may depend on 
the distance between a source and observation positions. To study position- 
dependent probabilities, the standard S-matrix of plane waves that satisfies 
a boundary condition at t — ±oo, which gives the values at the asymp- 
totic regions, is useless. An S-matrix that satisfies boundary conditions at 
a finite-time interval T, S[T], and has a position dependence is appropriate. 
Because boundary conditions for S[T] are different from those of 5[oo], it 
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reveals different properties form those of S"[oo]. We compute the finite-size 
corrections of transition probabilities with S[T] expressed by wave packets. 
The wave packets vanish at a position x if the distance between x and the 
center position X, \x — X\ — > oo, and satisfy the boundary conditions of the 
experiments [HE]. Hence they are appropriate to study the finite-size effect 
and are used here [TQ1 [HI [12] ■ 

A pion is produced in a proton reaction first and decays next. The whole 
process is expressed in Fig. 1 of I. The first reaction caused by strong inter- 
actions was studied in I and the second reaction caused by weak interactions 
is studied in II. From I, a particle in the final state retains the wave nature 
in the region of r < l , where r is a distance between the initial sates and 
final state and a coherence length / in a high-momentum region, \p\ 3> m, 
is given by 

<<=m- 

where m is the particle's mass and becomes for a neutrino of mass 1 [eV/c 2 ] 
and energy 1 [GeV] 

^eutrino = _± x K^eV" 1 ] = 10 2 - 10 3 [m], (2) 

Thus Zo °f neutrino are macroscopic lengths. Hence a measurement of neu- 
trino process at a near-detector region may depend on a distance \X — X^\, 
where X and X® are the positions of the nucleus in the detector and target. 

Detecting neutrinos at X of a distance \X — X^'\ Iq is studied in a 
usual manner. Since the neutrino in this region is expressed by a plane wave 
that behaves like a free particle, the neutrino's number and flux behave like 
those of classical particles and the production and detection of the neutrino 
are treated separately. The neutrino flux is determined by its distribution 
functions of the number and velocity that are determined by the decay pro- 
cess and is used for calculation of the scattering processes. 

Measurements in \X — X^'\ < Iq, on the other hand, are not treated by 
particle dynamics but by wave dynamics. The entire processes should be 
taken into account and be studied either by time-dependent Schrodinger 
equations or by position-dependent amplitudes defined with wave packets. 
The wave packets are necessary to measure the quantities in this region. 
Since the particles are identified based on signals generated by their interac- 
tions with a nucleus, atom, or larger system of matter, this unit of detector 
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has a finite size and is expressed by wave packets. Using wave packet repre- 
sentations, the amplitude of detecting neutrinos in a pion decay is computed3- 
The probability to detect the neutrino depends on the distance \X — X^\, 
hence a naive neutrino flux is not defined uniquely. The uniform flux defined 
by the classical particles is necessarily constant and is not used here. It might 
be reasonable to define the flux here based on the number of events that the 
neutrino gives rise to. The probabilities can be computed relatively easily 
with a set of simple wave packets such as Gaussian wave packets and are 
equivalent in more general wave packets. Although they depend upon not 
only the distance but also the wave packet size, they have universal proper- 
ties. So it would make sense to define a neutrino flux from the probability 
computed with the Gaussian wave packet. Thus we define the neutrino flux 
in this region from the probability to detect the neutrino. The idea is similar 
to define the electric field from a force that a point particle receives. The 
finite-size effects are actually large in the processes of detecting neutrinos 
at a macroscopic distance. Thus the corrections are important when the 
theoretical values are compared with the experimental values. 

It is instructive to see how a single neutrino interference is observed. A 
neutrino produced in a pion decay propagates a finite distance before it is 
detected. A position where a neutrino is produced varies, so a neutrino wave 
at the detector is a superposition of those waves that are produced at different 
space-time positions. If these space-time positions are inside of one pion, as 
in Fig.[TJ the waves keep their coherence and reveal interference patterns. 
A condition for the interference phenomenon to occur, for one dimensional 
motion of the pion that keeps a coherence within the region y / ov and a 
velocity v n is obtained in the following manner. Let a neutrino be produced 
either at time t\ or t 2 from the pion prepared at T^ and travel for some 
period and be finally detected at T„, then the waves overlap if 

\(c(T v - £x) + u w (tj - TV)) - (c(T„ - t 2 ) + v n (t 2 - TV))| < v^W, (3) 

is met, where the speed of light is used for the speed of neutrino v v = c. So 
Eq. (J3J) is one of the necessary conditions for the neutrino interference to occur 
in the one- dimensional space. For a plane wave of pion a n = oo and the above 
condition is satisfied. For a high-energy pion of a finite a n , its speed is close to 

1 The general arguments about the wave packet scattering are given in [13l 1X4} \TE \ I16 j . 
In these works, large wave packets were considered and small wave packets and finite-size 
corrections are studied here. 
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Fig. 1: The geometry of the pion decay region and the neutrino detector. 
The neutrino is produced in the pion decay and is detected. Since the decay 
position is not fixed to one value but is arbitrary, the amplitude to detect 
the neutrino is the overlap between the superposed initial wave and a final 
state expressed by the wave packet of the small size. The neutrino observed 
by the detector at T v shows an interference pattern. 



the speed of light and the left hand side of Eq. ([3]) becomes c{w?J r lE <2 ^){t\— t-i). 
Hence this condition Eq. ((3]) is written in the form, c{t\—t2) < ^fo\ ! {2E <2 K jrn l v ). 
When this length c(ti— £2) is a macroscopic size, the interference phenomenon 
occurs at a macroscopic length. We estimated the lengths of these particles in 
Appendix of I and confirmed that this condition in three-dimensional space 
is fulfilled in a macroscopic distance. 

An amplitude to detect a neutrino is an integral of a product of wave func- 
tions of a pion, muon, and neutrino and a velocity of a space-time position 
of transition of the pion to the muon has a component of the light velocity 
and the neutrino has almost the light velocity. Hence they overlap in wide 
space-time area and the amplitude and probability to detect the neutrino 
get a contribution of a large interference effect. A neutrino flux becomes 
distinct from the naive value obtained from the standard method of using an 
incoherent probability. 

This paper is organized in the following manner. In section 2, we study 
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an amplitude of detecting a neutrino in a pion decay process and compute 
a position-dependent probability in section 3. For a rigorous calculation 
of the position-dependent probability, a correlation function is introduced. 
Using an expression with the correlation function and its singular structure 
at a light-cone region, the finite-size correction is computed in section 4. 
Implication to neutrino experiments, features of the finite-size corrections, 
and summary and prospects are given in section 5, 6, and 7. 



2 Position-dependent amplitude of neutrino 

Position-dependent probabilities of detecting a neutrino and a charged lepton 
in a pion decay process are computed with a Schrodinger equation at a finite 
time t and an S-matrix of a finite-time interval T. They have various prop- 
erties different from those at the infinite time. Especially the final states in 
ultraviolet energy regions couple and give a finite and universal contribution. 

A neutrino or a charged lepton is detected with a detector that is located 
in a macroscopic distance from the position of the initial pion. Since their 
properties are determined from their probabilities in experiments, the whole 
process is studied. The process is described with a many-body wave function 
of the pion and decay products. The neutrino and charged lepton in the final 
state is highly correlated with the pion in the initial state and have a kinetic 
energy different from the initial value at a finite t due to a finite interaction 
energy. Hence the many-body state retains wave natures and a finite-size 
correction to the probability to detect a neutrino is expected. The scattering 
matrix of the finite-time interval S[T] is applied and a finite-size correction is 
computed using the wave packet representation. We find that the position- 
dependent probability is computed in a unique manner and reveals unusual 
finite-size correction which is understood as neutrino diffraction. 

2.1 Leptonic decay of the pion 

A leptonic decay of a pion is described with the weak Hamiltonian 

H 1 = g J dxd fl (p(x)J^_ A (x) = -igm^ J dxip(x)J 5 (x), 
J£_ A (x) = ft(x)^(l - lb )v{x), J 5 (x) = ft(x)(l - lb )v{x), (4) 
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where f(x), fi(x), and v{x) are the pion field, muon field, and neutrino field 
[T7| . In the above equations, the interaction Hamiltonian was expressed in 
a form of a derivative of the pion field, g is the coupling strength, Jy_ A (x), 
and J*,{x) are a leptonic charged V — A current, and a leptonic pseudoscalar. 
The coupling strength is expressed with Fermi coupling Gf and a pion decay 
constant f n , 

An expression of the weak Hamiltonian in term of (V — A) x (V — A) charged 
current interaction gives the equivalent results to a muon mode. In section 
6, an electron mode will be studied. 

A pion decay process is described by a Schrodinger equation 

ih^Mt)) = H\V(t)}, H = H + H l7 (6) 
where H Q stands for the free Hamiltonian. The solution is 

|*(t)> = e- 4 f*|*(0)> ) (7) 

and satisfies 

ih^wwmt)) = mm, mm)) = o, (s) 

ih^m)\H m)) = (m\iHo,H)m)) ? o. (9) 

The total energy is conserved but the kinetic energy is not conserved except 
the case 

H \V(t)) = E\V(t)). (10) 
A wave function of the initial condition 

Mt))\t=T„ = |pion(t))| t=T ^ = e-^lp^Tj, (11) 



evolves with 



|*(t)) = e - < f t |*(0)), (12) 
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and is studied hereafter. A solution in the first order in H% is a superposition 
of a pion state and a muon and neutrino state 

|\I/(t)) = |pion(t)) + |muon, neutrino(t)), (13) 

and each component is expressed as 

|pion(t)) = a(i)|pion, $,-(£)), (14) 

r f dt' 

|muon, neutrino(t)) = / — i?i(t')|pion(t')), (15) 

where a(t) = 1 + 0(g 2 ). The wave function Eq. ffl5|) is a superposition of 
states and is written as 



| muon, neutrino (t)) = ge 1 * l J dp^dp u ^ ^ M v 



e -iwt/h _ y 

X 5 (3) fe ~Vii -Vv) X (p*) li ji(Pr)'f(l ~ l5MPv)\P^Pv}> (16) 

where uj = + E u — E n and \pn,p v ) is a two-particle state composed of the 
muon and neutrino of momenta p^ and p u . 

At the infinite t, the oscillation function is approximated with 

p —iu)t/h i 

= -2iri8(u), (17) 

and the two particle state has the energy 

■ E I I 777 777, 

|muon, neutrino(t)) = -ige~ 1 ^ 1 / dp^dp v A——^—^-— (18) 

J V E n{Pl*) E »{Pv) 

x (2tt)(5 (4) (^ -p M -p„) x (^) At/ u(p M ) 7 ^(l - ^vip^p^Pu). 
The norm of this state is proportional to T and is given by, 

(muon, neutrino (T) | muon, neutrino(T)) = TT, (19) 



where V is the average decay rate [TU 19J. A neutrino and muon exist with 
a same probability that is computed from the norm of wave function. We 
should bear in mind that this is the probability that the decay products exist 
and a connection with an observation is not clear. Since an observation of 
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a neutrino is made with a detector that is located at a different position 
and a distance between positions of the pion and detector is usually large, 
the probability to detect the neutrino at the detector is defined according a 
geometry of the experiments with the wave function at a finite t, which is 
affected by a finite time and retarded effect. 

At a finite t, the oscillation function is broad in uj and Eq. ( II 7p does not 
hold 

-iuit/h _ i 

+ -2mS{u), (20) 

UJ 

and the two particle state Eq. (1161) has the continuous kinetic energy E^ + E u . 
The wave of the entire process in this region is described by 

| pion(t)) + |muon, neutrino(t)). (21) 

Because the two particle state composed of the muon and neutrino at a finite 
t has continuous kinetic energy, this state is different from two free particles, 
which has a constant kinetic energy, and is similar to the waves around a 
potential of finite energy. The wave retains a wave nature that varies with a 
time and the probability to detect a particle in the final state should reflect 
this wave nature. The S-matrix of observing the particle at a finite time 
t satisfies the boundary condition at a finite-time interval T, hence this is 
different from the standard S-matrix £[00]. We write this as S[T]. S[oo] 
commutes with H but S[T] satisfies, from the proof given in I, 

[S[T] , H ] = i ( ( t)) 1 Q + (T) - in! (T) (T ) . (22) 

Hence the states of non- conserving kinetic energy couple. The neutrino is 
measured from its collisions with nucleus in a detector, hence the wave func- 
tion of this final state in a scattering process is a wave packet of the nuclear 
size. Using the wave packets, we express S[T] following a reduction formula 
of [7]. This scattering amplitude depends on the position as well as the mo- 
mentum. The physics of neutrinos in this area is studied with the complete 
set of wave packets [10]. Since all the informations of the wave function 
I muon, neutrino (t)) are included in matrix elements, as far as a complete 
set of l/i, v) is used, these matrix elements have a complete information. 
Wave packets have been applied to study neutrino flavor oscillations before 
[2Ql EH [221 ESI EH [261 E7] , whereas S[T] and the probabilities that depend 
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on T are studied in the present paper. The finite-size correction appears in 
systems of one flavor and many flavors. 

Using a wave function of an initial pion located at a position X n , a neu- 
trino at a position X v and a muon, we express the amplitude to detect a 
neutrino at a finite distance. Since a muon is not observed and all the states 
are summed over, that is expressed with a plane wave. This transition am- 
plitude is 

T = J d A x{ii,v\H x {x)\<K), (23) 

where the pion and neutrino states are described in terms of wave packets of 
central values of momenta and coordinates and the widths in the form 

\k) = \p^,X^,T^), \n,v) = l^p^^p^X^T,,). (24) 

The pion is also expressed in term of a wave packet here for a completeness, 
and the particle states are defined with the matrix elements, 

(0| ¥ >(a;)|£ r ,.£ r ,T, r ) = N 7T p 7T J djg )re -¥(^-A) a - i ( s (^)( t - T -)-^-(^)) 

3 

(fi, v, p v , X v , T v \p(x)(l - 75) v(ac) |0) 



i 



A " ^Ke-^-^ 2 f i ' J£- ] ^)(1- T5 )^) 



where 



w '=(T)'-*=(T) ! ^=(^b)'- < 27 > 

In the above equation, the pion's life time is ignored. The sizes, a n and 
a v , in Eqs. f l25l) and ( |26i) are sizes of the pion wave packet and of the neutrino 
wave packet. Minimum wave packets are used in majorities of the present 
paper but non-minimum wave packets are studied and it is shown that main 
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results are the sameU From the result of Appendix of I, the coherence length 
of the pion produced in proton nucleon collision is short hence the pion is in 
the asymptotic region in a wide area. The pion behaves like a free particle 
there and can be treated with a plane wave or a wave packet. A mean free 
path of the pion was estimated in the Appendix of I and is used as a size of 
wave packet of the initial pion in a decay process in II. 

From Appendix of I, a size of pion wave packet is of the order of 0.5 — 
1.0 [m] and a momentum spreading is small. So k n is integrated easily, and 
is replaced with its central value p n and the final expression of Eq. (|25|) is 
obtained. For neutrinos, a size of wave packet is the nucleus. Hence to study 
neutrino interferences, we use the nuclear size for o v . 

The amplitude T for one pion to decay into a neutrino and a muon is 
written in the form 

T = J dtdxT(t,x), T(t,x) = J dk u T(t,x,k u ), (28) 

where 

T(t,x, k v ) = ^gr^A^'e"^^ - ^"^*"^) 2 "*^^* - ^^ 

x u(pj(l - 75 ) z/ (4)e i ( s ^)(*- T ^-^-^))-^^-^) 2 , 

3 1 



N^N.N^-J Ar„,*.= ^ x^- (29) 

T(t, x, ky) depends upon the coordinates (t, x) explicitly and is not invariant 
under the translation. 

The amplitude Eq. ( 12 8 p has the term from the region of k v ~ p„, which 
is computed with the S-matrix of plane waves, and another term from the 
region of k v ^ p u , in which the energy is not conserved. The latter term is not 
computed with the S-matrix of plane waves and is the finite-size correction 
that vanishes at T — > oo. We compute both with S[T] in the following. 



2 For non-minimal wave packets which have larger uncertainties, Hermite polynomials 
of k v — p„ are multiplied to the right-hand side of Eq. (l26l) . A completeness of the wave 
packet states is also satisfied for the non-minimum case |10) and a total probability and a 
probability at a finite distance and time are the same as far as the wave packet is almost 
symmetric. This condition is guaranteed in the high energy neutrino which this paper 
studies, but may not be so in the low energy neutrino. We will confirm in the text and 
appendix that the universal long-range correlation of the present work is independent of 
the wave packet shape as far as the wave packet is invariant under the time inversions. Low 
energy neutrinos such as solar or reactor neutrinos will be presented in the next paper. 
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2.2 Wave of observed neutrino : small angular velocity 
of a center motion 

Integrating the neutrino momentum in Eq. (|28p with the Gaussian integral, 
we find a coordinate representation of the neutrino wave function. 

For a not so large \t — T u \ region, k v is integrated around the central 
momentum p u , and the integrand T(t,x) becomes, 

x«fe)(l - 7 s)f (P^e^e-^^--^^"^, (30) 

where N is a normalization factor, v v is a neutrino velocity, and is a phase 
of neutrino wave function. They are given in the form 

N = N n Nj-Y (-Y N 0j (31) 



<t>{x) = E{p v ){t-T v )-p v -{x-X v ). (32) 

The neutrino wave function evolves with time in a specific manner. At 
t = T v , the wave is given in the form 

ij u (T»,x) = e^-akP-*") 3 , (33) 

which is localized around the position X v and has the phase = —p v -(x—X v ). 
At a time t < T v , the wave function becomes 

x) = e^-zk^-^-Mt-Tv)) 2 ^ ^34) 
which is localized around the position 

x G = X v + v v (t-T v ), (35) 

and has the phase 4>(x). Now this phase <p(x) is written at a position f = 
x — xq in the form, 

0(x) = G + 0(r), (36) 



where 

4>G 



- T„) - p„ • - T„) (37) 
3^AzlLa — r v\— m - n _ T ^ 

EM) 1 M ~*U^r " J ' 

— p ■ r. (38) 
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A phase at the center, <pG, has a typical form of the relativistic particle and 
is proportional to the mass squared and inversely proportional to the energy. 
Since the position is moving with the time, the phase from both components 
are cancelled and becomes extremely small. 

When the position is moving with the light velocity in the parallel direc- 
tion of the momentum p v instead of Eq. ( 13 5 p 



4> c (t - T„) = E{p v ){t - TV) - vw c(t - T„) = OF "A t - T„), (40) 



and becomes a half of (pe- 
lf the coordinate r is integrated in the amplitude, the phase 4>(f) is com- 
bined with those of the pion and muon fields and disappears at the end. 
We will see that in a process of detecting the neutrino with a detector at a 
long distance, an interference phenomenon of waves due to this slow phase 
occurs. This interference is almost equivalent to diffraction of light through 
a hole of a finite size. A light observed at a screen perpendicular to the wave 
vector shows a diffraction pattern, while a diffraction pattern is formed in a 
parallel direction to the momentum of the neutrino produced in the decay 
of pions. Since angular velocities of these phases are extremely slow and in- 
versely proportional to the neutrino energy, they make a diffraction pattern 
macroscopic. From a weak energy dependence at a high energy, the phase 
at the energy E v and E v + AE U is almost identical and is stable under a 
change of the initial state. So the neutrino diffraction has the same pattern 
in a broad energy spectrum and has the stability. 

The slow phase <p c shows a characteristic feature of the neutrino wave 
packet, and is an intrinsic property of the neutrino wave function at the light 
cone. This phenomenon is independent from the detail of wave packet and 
it is shown that wave packets in general cases, including those of spreading 
wave packets, have the phase. Since the phase in the longitudinal direction 
is the origin of the diffraction, spreading of the wave packet does not change 
the behavior in the longitudinal direction and does not modify the diffraction 
pattern. So the spreading effect has been ignored for simplicity in this section 
and will be studied in the latter section and Appendix. It will be shown there 




(39) 



the phase is given by 



ml 
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that the spreading in the transverse direction modifies the k u integration but 
the final result turns actually into the same. 

It is worthwhile to compare a neutrino velocity with the light velocity for 
a later convenience. A neutrino of energy 1 [GeV] and a mass 1 [eV/c 2 ] has 
a velocity 

v/c=l- 2e, e = (^jf-^ ' = 5 x 10" 19 , (41) 

hence the neutrino propagates a distance I, where 

I = l (l - e ) = l - 51, SI = el , (42) 

while the light propagates a distance Iq. This difference of distances, 81, 
becomes 

SI = 5 x 1(T 17 [m]; l = 100 [m], (43) 

SI = 5 x 10~ 16 [m]; l Q = 1000 [m], (44) 

which are much smaller than the sizes of the neutrino wave packets. Since 
the difference of velocities is small, the neutrino amplitude at the nuclear or 
atom targets show interference. The geometry of the neutrino interference 
is shown in Fig. [TJ The neutrino wave produced at a time t\ arrives at one 
nucleus or atom in the detector and is superposed to the wave produced 
at t 2 and arrives to the same nucleus or atom same time. A constructive 
interference of waves is shown in the text. 



3 Position-dependent probability and inter- 
ference 

Because the wave function at a finite time is a superposition of waves of 
broad kinetic energies, that retains the wave natures and reveals a diffraction 
phenomenon. The probability to detect the neutrino reflects this property 
and receives a large finite-size correction peculiar to tiny masses in a wide 
area. The finite-size correction is computed rigorously with a use of the light- 
cone singularity of a two point correlation function of a decay process. The 
correction has a universal property of the relativistically invariant system 
and is determined by the absolute neutrino mass. 
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3.1 Transition probability 

We investigate the case of large ov later and study the amplitude and prob- 
ability of ov = oo qualitatively first. Features of the amplitude to detect the 
neutrino at a finite distance are elucidated. Integrating over the momentum 

— # 

k in Eq. (|28p . we obtain the ^-dependent amplitude of the form Eq. (15U|) . The 
integrand is a Gaussian function around the center x (t) = v u {t — T u ) + x(0) 
and is invariant under 

x ->■ x + v v 5t, (45) 
t->t + 6t. 

Thus a shifted energy 

H -v w - P, (46) 

satisfies 

[S,H -v u -P}=0 (47) 

and is conserved. 

Integrating over x further, we have the amplitude, 

T = Ce^uip.)(l - j 5 )u{p v )e-^ 2 e-^ T / 2 2 ^ , (48) 

to 

u = 5E -v v - 6p, 8p = p v -pp- Pu, SE = E{p v ) - E(p /J ) - E{p v ), 

where 0o is a constant. Because the center x*o{t) moves with the velocity v u , 
the angular velocity in Eq. ( 148 p is different from the energy difference 5E of 
the rest system, but 8E —v v -5p of the moving system. This is a feature of 
the present amplitude. 

In Eq. (jUJ), the momentum is approximately conserved due to the Gaus- 
sian factor and \8p \ has a finite uncertainty. Hence the angular velocity u> 
behaves differently from SE. At Sp = 0, u = is the same as the usual case 
SE = 0, whereas at Sp ^ 0, uj = gives the relation SE — v v -Sp ^ 0. Kinetic 
energy takes broad range because u is different from SE and the amplitude 
is broad in u at a finite T. Thus, kinetic energy is not conserved SE ^ 0. 
As is shown in Appendix B, a shape of the configuration of the momentum 
satisfying to = is a large ellipse of the muon momentum where the nor- 
mal solution of Sp = and the solution of large \Sp\ are on the curve, u 
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varies fast around the former momentum and 2 S1D ^ T ^ 2 ' > = 2tt5(uj) [28] can be 
applied. This gives the normal term which satisfies the energy-momentum 
conservation well. On the other hand, u varies extremely slowly around the 
latter momentum, and the states of u ~ lead the slow convergence at large 
T, and give the finite-size correction. Since \5p\ and 5E are not small, the 
spectrum at the ultraviolet region, which exists in the wave function at a 
finite time, gives a contribution to the finite-size correction. We will study 
this point further in Appendix B. 

Although the reason for the large finite-size correction became clear, it 
is not straightforward to compute it using Eq. (148]) . Instead, an expression 
of |T| 2 with a correlation function of coordinates is convenient for a rigorous 
computation, and the finite-size correction is found with it. 

A transition probability of a pion of a momentum p^ located at a space- 
time position (TV, X n ), decaying to the neutrino of the momentum p v at a 
space-time position (T u , X u ) and a muon of momentum p M , is expressed in 
the form 



|T| 2 = g 2 ml\N\ 2 J d 4 x l d 4 x 2 S 5 (s u s 2 ) 



x e -i(E(pTr)(t 1 -t 2 )-p 7T -{x 1 -x 2 )) x e i{E{p li )(t 1 -t 2 )-p^{x 1 -x 2 )) 

where S^Si, s 2 ) stands for products of Dirac spinors and their complex con- 
jugates, 

S 5 (si, s 2 ) = («(p At )(l - 75)f(p«/)) («(Pm)( 1 ~ 75)K&))* > ( 50 ) 
and its spin summation is given by 

^ 5 = ^5 5 ( Sl , S2 ) = ^— (p,- Pu ). (51) 

The probability is finite and an order of integrations are interchangeable. 
Integrating momenta of the final state and taking average over the initial 
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momentum, we have the total probability in the form 



dpnPexpiPVjdXudpudpv ^2 \ T \ 



Nm = O 4 © 4 ' pv = (^f) 2 ' 6x = xi - x2 > (52) 

with a correlation function A ntfl (St,Sx). The correlation function is defined 
with a pion's momentum distribution p exp (p n ), by 

(53) 

In the above equation, the final states are integrated over a complete set 
[TO] . The muon and neutrino momenta are integrated over entire positive 
energy regions, and the neutrino position is integrated over the region of 
the detector. The pion in the initial state is assumed to be the statistical 
ensemble of the distribution p e xp{Pn)- If the momentum distribution is narrow 
around the central value, the velocity v n in the pion Gaussian factor was 
replaced with its average v n . This is verified from the large spatial size of 
the pion wave packet discussed in the previous section. For the probability 
of a fixed pion momentum, the correlation function 

^ ' ° X) ~ (2tt)3 Efa) J E(p„) ^ Vv) 

(54) 



is used instead of Eq. ([53 



3.2 Light-cone singularity 

The expression Eq. (1521) shows that the probability gets a finite T(T = 
T u — T n ), correction from the integration over t\ and t 2 at \ti — t 2 \ — > T, 
if the A nttl (St, 8x) decreases slowly in this region. The correlation function 
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A n<ll (St,Sx) is a standard form of Green's function and has the light-cone 
singularity that is real and decreases very slowly along the light cone. The 
singularity is generated by the states at the ultraviolet energy region. So the 
singularity near the light-cone region 

A = 5t 2 - \5x\ 2 = 0, (55) 

which is extended in a large \5x\ and is independent of p n , plays an important 
role for the probability Eq. f )52|) at a finite T. We find, in fact, that the light- 
cone singularity of A x ^(5t, 5x) gives a large finite-size correction in the 
following. 



3.2.1 Separation of singularity 

If the particles are expressed by plane waves, the integration is made over 
infinite-time interval and the energy is strictly conserved and 4-dimensional 
momenta satisfy 

fV = Vii + Pu, (Pn ~ P^Y = m 2 u ^0. (56) 

An integral over the momentum in the region where the momentum difference 
Pw—Pfi is almost light-like leads A n ^(5t, Sx) to have a singularity around the 
light cone, A = 0. In order to extract the singular term from A n tfi (8t, 5x), we 
write the integral in a four- dimensional form 

A^(5t,8x) = 1 * I(p T ,5x), 
(2ny E(p 7T ) 



I(p n ,5x) = -^f d 4 Pn°(pl)(Pn-P») lm 



(57) 



first, and change the integration variable from p^ to q = p^ — p n that is 
conjugate to 5x. Next, we separate the integration region into two parts, 
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< q° and — p° < q° < 0, and have the expressions, 



I(p w , Sx) = h(p w , Sx) + h(Pir, Sx) 

. d 



(58) 



dSx 
1 



h(pn,5x) -- 

h = - f A%°)Im , 

2 r° 

I 2 (p*, Sx) = - / (Pqp v -(p n + g)Im 



m,, — ie 



D iq-5x 



mf, — 



Jq-Sx 



(59) 



A(Ptt5^) is the integral over the infinite region and has the light-cone sin- 
gularity and ^{Pn, Sx) is the integral over the finite region and is regular. 

h{Pn, Sx) comes from the states of non-conserving kinetic energy and does 
not contribute to the total probability at an infinite-time interval. ^(Ptt, Sx), 
on the other hand, contributes to that at the infinite-time and finite-time 
intervals. So the physical quantity at the finite distance is computed using 
the most singular term of I\. 

Next we compute I\. Expanding the integrand with p^-q, we have I\ in 
the form 



<fq6{q 



2 

7T 

4„ /i/„0 



A%°) Im 



q 2 + m 2 — m 2 + 2g-p„ 



K 



D iq-6x 



= 2 
where 



2jv 



d 



2Ptt- 












dSx J dm 



Im 



g 2 + m 2 



IE 



iq-Sx 



dSx J dm 2 



+ 



d 4 q6(q°)S(q 2 + m 2 )e iq - 5x 



m 2 



ml - ml. 



The expansion in 2p 7T ■ q of Eq. (160]) converges in the region 

2p w -q 



q z + m 



- , 1 2 



< l. 



(60) 



(61) 



(62) 



Here q is the integration variable and varies. So we evaluate the series after 
the integral and find a condition for its convergence. We find later that the 
series after the momentum integration converges in the region 2p z'£ v < 1. 
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hiPm $ x ) is written in the form 



h(pn,6x) = 2(2tt) 3 ^ l + 2p w - [i 



d 



d 



d5x J dm 2 



+ 



47T 



S(X)e(St) + f short j , 
(63) 



where f short is written by Bessel functions and a formula for a relativistic 
field 



/ 



d 4 q0(q°)5(q 2 +rh 2 )e iq - 5x = (2n) 3 t 
im 



—5(X)e(5t) + f short 



short 



im 



0(-A) jiVi (mV^x) - ie(5t) Ji (mv^x) } 



- 0(A) ^^K x (rhVx) , X = 5x 2 , 5t = 5x°, 



■WT" "" -' (64> 

where Ni, Ji, and are Bessel functions, is used. More details are presented 
in I. 

Next I 2 is evaluated. For I 2 , we use a momentum q = q + p n and write 
in the form 



I 2 (p n ,5x) = - d 4 q(p u -q)lm 

n J0<q°<p° 

. d M 2 



g 2 — m 2 — «e 



3 i(g-p w )-6a; 



e~ ipiT ' 5x <,p v -\ -%- 



d5x 
d 



TV 



d 4 q~7i5(q 2 -mly 



0<g°<p£ 



(pl - sjq 2 + ml) e" 



iq-8x 



(65) 



85x J J J ^Jq 2 + m 2 

The regular part /2 has no singularity because the integration domain is finite 
and becomes short-range. 

Thus the first term in I\ gives the most singular term and the rests, 
the second term in l\ and I 2 , give regular terms. The correlation function, 
A,r jM ((ft, 5x) is written in the form 



A nilx (St,Sx) = 



d_ 

d5x 



2(2vr) 3 i 



l + 2p n -(i 



d 



d 



+ 



dSx J dm 2 

where the dots stand for the higher order terms 



— 5(X)e(5t) + f short ) + h 



(66) 
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3.3 Integration over spatial coordinates 

Next, we integrate over the coordinates x\ and xi in 

J rffxrffae^e-^^^-^-^^-^^'A^^e). (67) 
3.3.1 Singular terms: long-range correlation 

The most singular term of A 7r /i (5t, Sx) is substituted, then Eq. f l6T|) becomes 
J 5{X ) = J rfx 1 df 2 e^ (fa) e-^ E ^^-^~^ {tl - T " ) ) 2 ^5(A)e(^), (68) 

and is computed easily using a center coordinate i? M = 1 2 2 and a relative 
coordinate r = x\ — x 2 - After the center coordinate R is integrated, Js(\) 
becomes the integral of the transverse and longitudinal component (rj>, rf) of 
the relative coordinates, 

e(5t)(a v n)* / df T dn e l ^-^ {? THn-v v stf)^ 5t 2 _ f 2 _ ^ (6Q) 

J 47T 

The transverse coordinate fr is integrated using the Dirac delta function and 
r\ is integrated next. Finally we have 

/ '±O v \ ... iiJSt) ^^St 2 

J S {\) = ( a ^ 2 Y]5t\ e ^ e 

« K^y^^)^- (TO) 

The next term of A 7r A1 (5t, of the form j-, in Eq. ( 1671) leads 

J 1/A = /'rff 1 rff 2 e^)e-^^( a;i -^-^^- T ^) 2 ^-, (71) 
y 47r i A 



which becomes 



Ji/x « Ktt)^ — ,2 e"-^l TTire^. (72) 

2 \7T(7 I/ |p I/ |' 4 / |d£| 

This term also has the universal \St\ dependence but its magnitude is much 
smaller than that of Jjm and is negligible in the present decay mode. 
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From Eqs. (170]) and ( T72|) . the singular terms J^a) and J\/\ have the slow 
phase c (5t) and the magnitudes that are inversely proportional to <5t. Thus 
these terms are long-range with the small angular velocity and are insensitive 
to the ifi 2 . These properties of the time-dependent correlation functions JsiX) 
hold for the general wave packets, and the following theorem is proved. 

Theorem 

The singular part J^m of the correlation function has the slow phase that 
is determined with the absolute value of the neutrino mass and the magnitude 
inversely proportional to 5t, of the form Eq. (!70|) . at the large distance. The 
phase is given in the form of a sum of <f) c {5t) and small corrections, which are 
inversely proportional to the neutrino energy in general systems and become 
1/E 2 if the neutrino wave packet is invariant under the time inversion. 
(Proof: General cases including spreading of wave packet ) 
We prove the theorem for general wave packets. Js(\) is written in the 
form, 

J S(X ) = J dre l ^ Sx) w{r-v u 5t) X -jU(A)e(<ft), (73) 

where w(x — vt) is expressed with a wave packet in the coordinate represen- 
tation w(x — vt) and its complex conjugate as, 

w( ri - vjt, f T ) = j dRw (r + w* - 

d^^re^ (n ~^ 5t)+ ^ T ' f ' T+ico( ^)' 5 *|u;(A; z , k T )\ 2 . (74) 

The wave function w(x — vt) that includes the spreading effect is expressed 
in the following form 



w ( 



(f - m) = / dkidhr e ik ^- v ^ +i ^- ST + iCi ^T k ^w(k h k T ), (75) 



Cij - C 5ij, C — — , (76) 

instead of the Gaussian function of Eq. (168 p . A quadratic form in k in an 
expansion of E (p + itj is included and this makes the wave packet spread 
with time. The coefficient CV, in the longitudinal direction is negligible for 
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the neutrino and is neglected. Expanding the delta function in the form, 

5(5t* - rf - f|) = l(-4Y (^) ' *(* 2 - rf ), (77) 
we have the correlation function 

■Mx) = J dndrre^^win - v v 8t,r T )^ j 1 + E ^~^T ( & 
x 5(5t 2 - rf)e{St) 



d(sty 



T)\ 2 



x i I 1 + E ^H?)" {JL ) \ W> - rfHSt) 



4ir '—'nV " \d(St) 

\ n=l \ \ / 

= J dndhe tm ' rMikl{ri ~ VvH) "dr T dk T e +iCoP T St \w(k h k 



/Ml 2 



(78) 



The variable r T is integrated first and £>r is integrated next. Then we have 
the expression 

Js(\)= J rfn^e^^^^^-^^lwl^.o)! 2 

x I {l + E ^(-«WT (^) "} (2^)^(« 2 - -fM*). (79) 



Using the following identity 



\ I a \ „ / 1 \ / a 



i2Sty '{mm) ={m) +0 {T t ){m) • (80) 

and taking a leading term in we have the final expression of the corre- 
lation function at the long-distance region 



: 4c(St) 



Js(x } = ^ e ~ C ° P < St )^f- J dhe^'^^wihM 2 - (81) 
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Hence Jg(X) m Eq. ( 1ST]) becomes almost the same form as Eq. (ITU!) and 
the slow phase <f> c (5t) is modified slightly and the magnitude that is inversely 
proportional to the time difference. J$(\) has the universal form for the 
general wave packets. By expanding the exponential factor and taking the 
quadratic term of the exponent, the above integral is written in the form 



/ 



dki{l + k(i(l - v„)6t + Co) + ^(ki(i(l - v v )5t + C )) 2 )\w(k h 0)| 2 
w (l + C Q d x + ^C 2 - (1 - v v fsA + i{d x {\ - v v )St + d 2 C Q {l - v u )6t), 



(82) 



where 



d\ d 2 1 di d 2 ( 1 \ 2 , \2z+2 

d 1 = — f dhh\w{k h 0)\ 2 1 d 2 = — [ dkikf\w(k h 0)\ 2 . (83) 
Wo J w J 

We substitute this expression into the correlation function and have 

i$ c {6t)(l+S) f 

J 6{ x) = ^e- c ^ujo{l + i)e{5t) — , w = I dk\w{k h 0)\ 2 , (84) 

In wave packets of time reversal invariance, \w(ki, 0)| 2 is the even function 
of k\. Hence d\ vanishes 

di = 0, (85) 

and the correction are 

S = 7 = - x (— V - (1 - v v )Ht\ (86) 

2 E 2 ' 2! \2E J v ; v ; 

Q.E.D. 

The light-cone region 5t 2 — \5x\ 2 = is so close to neutrino orbits that 
it gives a finite contribution to the integral Eq. fl67|) . Since the light-cone 
singularity is real, the integral is sensitive only to the slow neutrino phase 
and shows interference of the neutrino. This theorem is applied to quite 
general systems, where the neutrino interact with a nucleus in a target. 
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3.3.2 Regular terms: short-range correlation 

Next, we study regular terms of A ntfl (5t,5x) in Eq. (j67]l . Regular terms are 
short-range and the spreading effect is ignored and the Gaussian wave packet 
is studied. First term is / short m h and is composed of Bessel functions. We 
have 

Li is evaluated at a large \St\ in the form 

L x = (7ra u )h iE " 5t [ dre-^-^~^ 2 f short , r = x x - x 2 . 



Here the integration is made in the space-like region A < 0. We write 

ri =v„6t + r h (89) 

and rewrite A in the form 

A = 5t 2 - rf - f| = 5t 2 - (v„6t + hf - Tt ~ -2v v f t St - f] - f|. (90) 

The Li for the large \5t\ is written with these variables. Using the asymptotic 
expression of the Bessel functions, we have 

J 4tt 2 \2mJ 

3 

4 , 



3 

/ 1 \ 

X 



2v v fi\St\ + rf + r| 

The Gaussian integration around r T = 0, fj = — 2%(j v \p v \ give the asymptotic 
expression of L x at a large \&t\ 



T = r i(Eu-Wv\Vv)6t -a v \y v ? 1771 ( 77 \ 2 ( ) * p im-v/2iv CT ,,|pi,||<5t| 

1 4vr 2 v2mJ \Av u a v \p v \\5t\ ) 

(92) 

Obviously Li oscillates fast as e imc ^ 5t ^ where c\ is determined by \p v \ and 
a v and is short-range. The integration carried out with a different stationary 
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value of r\ which takes into account the last term in the right-hand side gives 

almost equivalent result. The integration in the time- like region, A > 0, 

i 

is carried in a similar manner and I\ decreases with time as e ~ mci l' 5 'l 3 and 
final result is almost the same as that of the space-like region. It is noted 
that the long-range term which appeared from the isolated 1/A singularity 
in Eq. (172 j) does not exist in Li in fact. The reason for its absence is that 
the Bessel function decreases much faster in the space- like region than 1/A 
and oscillates much faster than 1/A in the time-like region. Hence the long- 
range correlation is not generated from the I\ and the light-cone singularity 
S(X)e(6t) and 1/A are the only source of the long-range correlation. 
Second term of Eq. fl67|) is from J 2 , Eq. f )65|) . We have this term, L 2 , 

■3, 3 | ~ 

^(47rer I ,)2 7 — —L 2 , 



x e 



2p v -(p 1T -p v )(-Ka lJ 
dq 

— = ( 

2 x Jq 2 +m* 

-ov(i>-K—q—Vv) 2 Q ( £ 



(2vrr 

- y/ V 2 +mf l -v lJ -(p 7v -q-pv))St 



q 2 + 777^ 



(93) 



The angular velocity of the integrand in L 2 varies with q and the integral 
Z/2 has a short-range correlation of the length, 2^/d^, in the time direction. 
So the L 2 's contribution to the total probability comes from the small \5t\ 
region and corresponds to the short-range component. 

Thus the integral over the coordinates is written in the form 



dx x dx 2 e^ (fa) 



2i 



e -^E,^-i-4(t«-T,)) ^ 

d 



7T,/i. 



(St, Sx) 



1 + 2p n -p 



2i(na l/ 



dm 2 
1 + 2p 7T -p 



+ 



J4>{St) i 



{J s m + In) + U 



d 



x 



\5t\ 



+ Ii 



dm 2 



+ 



(94) 



In the above equation, p v = m v is negligibly small compared to m , p v -p v 
and (Tj, -1 , and is neglected in most places except the slow phase 4>(5t). The 
first term in the right-hand side of Eq. ( 194"|) is long-range and the second 
term is short-range. The long-range term is separated from others in a clear 
manner. 
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3.3.3 Convergence condition 

At the end of this section, we find a condition for our method to be valid. In 
Eq. f )58|) . the integration region was split into the one of finite region — p° < 
q° < and the region < q°. Accordingly, the correlation function is written 
into a sum of the singular term and the regular term. The singular term is 
written with the light-cone singularity and the power series in Eq. flBTJ]) . Hence 
this series must converge for the present method of extracting the light-cone 
singularity to be applicable. 
We study the power series 

using the asymptotic expression of Li, Eq. ( |92|) . first. The most weakly con- 
verging term in L 1; is from ma and other terms converge when this converges. 
The series 



m 2 )2. (96) 



becomes the form, 



-2p 7T -p u \ 1 i 1 , _ 
n ! ( — 1 (m 2 



m 2 I n! V i 



- E (-^)" (-!)•»-»(*)* = E (^)""- f ('»)'- < 97 ' 

Hence the series converges if the geometric ration is less than 1 . At 2^ • p u = 
rh 2 Si becomes finite, and the value is expressed by the zeta function, 



(98) 



Hence in the region, 



m 2 



< 1, (99) 



the series converges and correlation function has the singular terms. Outside 
this region, the power series diverges and our method does not work. / is 
evaluated directly and agree with the J 2 . 
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The power series Eq. ()95p oscillates with time \/\8t\ rapidly. The series 



1 / F) \ n 

s 2 = Y,h2 P7T -p u ) n -(-^) e^V^^m, (100) 

n ■ \ / 

is computed in the form, 

s 2 = g^lV^^i^n^iK 1 - 2 ^) (101) 

and oscillates with So the present method of separating the light- 

cone singularity from the correlation function and of evaluating the finite-size 
correction of the probability is valid in the kinematical region Eq. (199 p . 



3.4 Time-dependent probability 

Substituting Eq. ( 19^1 into Eq. (152]) . we have the probability for detecting 
the neutrino at a space-time position (T U ,X U ), when the pion momentum 
distribution p e xp{Pw) is known, in the following form 




(102) 

From a pion mean free path obtained in the Appendix of I, the coherence 
condition, Eq. ([3]), is satisfied and the pion Gaussian parts are regarded as 
constant in ti and t 2 , 

e -&(&-&+fr^m-T„)-&(?«-T u )) a « constant in t u 
e -^(^-^+(tt-^)(t 3 -T„)+^cr w -T„)) a w congtant in (1Q3) 

when an integration over t\ and t 2 are made in a distance of our interest 
which is of the order of a few 100 [m]. The integration over t\ and t 2 will be 
made in the next section. 
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When the above conditions Eq. (11031) are fulfilled, an area where the neu- 
trino is produced is inside of a same pion and neutrino waves are treated 
coherently and are capable of showing interference. In a much larger dis- 
tance where this condition is not satisfied, two positions can not be in the 
same pion and the interference disappears. 

3.4.1 Integrations over times 

Integrations over the times t± and t 2 are carried and a probability at a finite 
T is obtained here. An integral of the slowly decreasing term over t\ and t 2 
is 



/ dt 1 dt 2 —-e(8t) = T{~g(T,co v )+g(oo,Lo u )}, u u = — ^,(104) 
Jo \ot\ zh v 



where g(T,cov) vanishes at T — > oo. We understand that the short-range 
part Li cancels with g(oo,u u ) and write the total probability with g(T,u u ) 
and the short-range term from I 2 . 

An integral over times of the short-range term, L 2 , is 



7T V vr 



dtidt 2 L 2 (5t) 



1 K [ T dhdt 2 f , j5 „ e -<^-^-v^-g*-(ft-7-fto)« 

7T V 7T 







x e-^P^-P^e (e„ - yjq 2 + rr. 
= TG , (105) 

where the constant Go is given in the integral 

G o = 2^ J ^F=p 5 { E « ~ E « -\J^ + m l-^ ■ & - 9-Pu)) 

and is estimated numerically. Due to the rapid oscillation in 5t, L 2 s contribu- 
tion to the probability comes from the small \St\ region and the integrations 
over the time becomes constant in T. Hence this has no finite-size correction. 
The regular term Li is also the same. 
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Fig. 2: The probability to detect the neutrino per time in the forward di- 
rection at a distance L is given. The constant shows the short-range normal 
term and the long-range diffraction term is written on top of the normal 
term. The horizontal axis shows the distance in [m] and the probability of 
the normal term is normalized to 2tc. Clear excess of more than 2/5 of the 
normal term is seen in the distance below 1200 [m]. The neutrino mass, pion 
energy, neutrino energy are 1 [eV/c 2 ], 4 [GeV], and 800 [MeV]. Target is 16 0. 

3.4.2 Total transition probability 

Adding the slowly decreasing part and the short-range part, we have the 
final expression of the total probability. The neutrino coordinate X u is inte- 
grated in Eq. f 1 1 2 j) and a factor (a n 7r) a emerges. This factor is cancelled with 

3 

(47r/cr 7r )2 of the normalization in Eq. (149 ft and a final result is independent 
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of a n . The total transition probability is expressed in the form, 



P = Tg m^D a v J —p, 




1 



dpi 
E v 
1 



D = |A^|V^) 2 (°V7r) 2 



3 1 



(107) 



_> 



(27T) 



(2tt) 3 ' 



where L = cT is the length of decay region. The first term in the right- 
hand side of Eq. (11071) depends on the time interval T, and the neutrino wave 
packet size a u , but the second term does not. 

At a finite T, the first term, which we call a diffraction term, does not 
vanish and the probability Eq. f ll07f) has the finite-size correction. Its relative 
ratio over the normal term Go is independent of detection process. So we 
compute g(T, uj v ) and G of Eq. fllOTf) at the forward direction 9 = and the 
energy dependent total probability that is integrated over the neutrino angle 
in the following. 

The probabilities per unit time in the forward direction are plotted in 
Fig. [2] for the mass of neutrino, m„ = 1 [eV/c 2 ], the pion of the sharp energy 
E n = 4 [GeV], and the neutrino energy E v = 800 [MeV]. For the wave packet 
size of the neutrino, the size of the nucleus of the mass number A, a v = 
Az/rrin is used. The value becomes a v = 6.4/m 2 for the 16 nucleus and 
this is used for the following evaluations. From this figure it is seen that 
there is an excess of the flux at short distance region L < 600 [m] and the 
maximal excess is about 0.4 at L = 0. The slope at L = is determined by 
u u . The slowly decreasing term has the finite magnitude and the finite-size 
correction is large. 

4 Neutrino spectrum 

4.1 Integration over neutrino angle 

In Eq. p07p . the diffraction term g(T,u u ) has a different dependence on the 
angle from that of the normal term Go- In the normal term Go, the cosine of 
neutrino angle 9 is determined approximately from a mass-shell condition, 



because the energy and momentum conservation is approximately well sat- 
isfied in the normal term. Hence the product of the momenta is expressed 



(108) 
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with the masses 



2 2 

mi — m,, 

P*-P*= 2 ; (109) 
and the cosine of the angle satisfies 

m 2 — m 2 m 2 
l-cosg= JL -J!^-. (110) 

2|£r||&| 2|k| 2 v ; 

The cos 9 is very close to 1 in a high energy region. On the other hand, the 
diffraction component, g(T,uj u ) of Eq. (11071) . is present in the domain of the 
momenta Eq. (1991) i.e., in the kinematical region, 

2 2 
ffi — Tfl 

\p u \(E n - \&\) < p n -p u < *- 2 (Ill) 

Since the angular region of Eq. (II 111) is slightly different from Eq. (11091) and it 
is impossible to distinguish the latter from the former region experimentally, 
the neutrino angle is integrated. We integrate over the neutrino angle of both 
terms separately. We have the normal term, Go, in the form 
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K V 11 \ a v) J ^Jq 2 + m^ 

X5(E*-E V - ^/g 2 + m 2 ) 5 (3) fa-fr - $6(E V - y^ + m 2 ) 

dE„, (112) 



(2-)'"' [ml - m 2 \ i r ^-<»"-< 



\P*\ 



where 



2 2 2 2 

mi — m„ mi — mi 

2(^ + 1^1) 2{E n -\p lv \) 

and the Gaussian function is approximated by the delta function for the 
computational convenience. The angle is determined uniquely. 

We compute the diffraction term next. There are two cases depending on 
the minimum angle of satisfying the convergence condition Eq. (l9~9~|) . In the 
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first energy region, 

E,E„- |K - mj) 

\p*\\pA ' [ ' 

the convergence condition is satisfied in a limited region of the angle. We 
integrate over this region of the angle and have the diffraction term in the 
form 

= ^ J f dcosOjE^ - l^l^l cos^T,^) 

\Ptt\\Pv\ 

= 2n Je ^{4 W-^-^-IAII^jsfCT.o;,,). (115) 

Here the angle is very close to the former value but is not uniquely de- 
termined. In the second region, the convergence condition is satisfied in 
arbitrary angle, 

E W E U - Uml - m 2 u ) / N 

— 2 , \ ^ < -1, (116) 

\pA\pA 

and we have the diffraction term in the form 

^r{p^-Pu)g{^^ v ) 

f \Pu\ 2 d\pA f 1 
2n J — ^ — — J dcos9(E n E u - \p n \\p v \ cos9)g(T,u u ) 

4tt / dE u E w E 2 J(T,u u ), (117) 







In the diffraction term, the angle of the neutrino is different from that 
of the normal term. The angle dependences of the energy of normal and 
diffraction terms are given in Fig. 3. The angle is fixed to one value in the 
normal term and is in a continuous range in the diffraction term. 
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5 10 15 20 25 30 35 40 

Angle [degree] 

Fig. 3: The relation between the neutrino angle and energy is shown. The 
energy is determined uniquely with a value on the line in the normal com- 
ponent and takes a finite range under the line in the diffraction term. The 
energy of the pion is E v = 2 [GeV]. 



Finally we have the energy dependent probability 



dP 
dE„ 



Tg 2 mlD \ -^rP exv (Vv 



2vr 

\Pn\ 



ir (ml - ml) 



+ y (V, " E v>min ) ji (ml - ml) 2 - (E K E V - \Mp»\) 2 } 
+ e(K, m in - E v )2E k eAq(T, ou u ) 



118) 
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4.2 Neutrino spectrum 
4.2.1 Sharp pion momentum 

When the initial pion has a discrete momentum P^, the PexpiPn) is given as 

Pexpijh) = - P -*)i ( 119 ) 

and the probability is expressed in the form, 



dP „ o o„ 1 2tt 



dE u y " °K|PJ 



7r(m 2 - mj) 



0"„ 



+ ^ - £„, min ) - m 2 ) 2 - - \P«\\p v \)' 

-\-0 (Ev,min 

- E v )2E w El) g(T, u v 



(120) 



Eq. (11201) is independent of the position X n and an average over X w is 
easily made. The result is obviously the same as Eq. (11201) . 

The probability depends upon the momenta and the time interval T = 
T„ — T w . At T — > oo, g(T,u u ) vanishes and the probability per unit time, a 
decay rate, of the energy E n is given from the first term of Eq. (I120p as 



i 2tt , „ r K 



P/T = g* m lD —- w -Tr(ml - ml) / dE, 



U 



■> / 2 \ 2 

,2,^.2 TV / i M 



7T 



9'<r*tl l -Ji- (121) 



This value is independent of the wave packet size and is consistent with [24|. 
This, furthermore, agrees with the standard value obtained using the plane 
waves. Hence the large T limit of our result is equivalent to the known result 
obtained with the standard S-matrix. 



4.2.2 Position dependence 

The position dependence of the total probability per time P/T, Eq. (11181) 
is presented next. g(T,u) u ) varies with the distance L defined by L = cT, 
whereas G is constant. The probability depends upon the magnitude u u T 
and for the neutrino mass m u = 1.0[eV/c 2 ] and the pion energy 4 [GeV] 
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Fig. 4: The total probability per time integrated over the neutrino angle at 
a finite distance L is given. The constant shows the normal term and the 
diffraction term is written on top of the normal term. The horizontal axis 
shows the distance in [m] and the total probability is normalized to a unity at 
L = 0. The excess becomes less clear than the forward direction, but is seen 
in the distance below 1200 [m]. The neutrino mass, pion energy, neutrino 
energy are 1.0 [eV/c 2 ], 4 [GeV] and 4.5 [GeV] , and 800 [MeV]. Target is 16 0. 



and 4.5 [GeV] are given in Fig.HJ and for the smaller neutrino mass m u = 
0.6 [eV/c 2 ] is given in Fig.[5j g(T, u v ) of a lighter mass decreases more slowly 
with the distance than that of m v = 1 [eV/c 2 ]. A longer distance is necessary 
to see a signal if the neutrino mass is even smaller. For the detection of the 
muon neutrino, the neutrino energy should be larger than the muon mass, 
hence the experiment of the energy lower than 100 [MeV] is impossible. For 
this energy, the electron neutrino is used then. We present the total proba- 
bility for the lower energies next. The probability for m u = 1.0 [eV/c 2 ] with 
the energy 100 [MeV] is given in Fig.[6j The slowly decreasing component 
of the probability becomes more prominent with lower values. Hence to ob- 
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Fig. 5: The total probability integrated over the neutrino angle per time at 
a finite distance L is given. The constant shows the normal term and the 
diffraction term is written on top of the normal term. The horizontal axis 
shows the distance in [m] and the probability of the normal term is normalized 
to 0.8. Clear uniform excess is seen in the distance below 1200 [m]. The 
neutrino mass, pion energy, neutrino energy are 0.6 [eV/c 2 ], 4 [GeV], and 
800 [MeV]. Target is 16 0. 



serve this component, the experiment of the lower neutrino energy is more 
convenient. 

From Eq. (11071) . and g(T,u) u ) = at a large T, the typical length Iq of 
the diffraction term is 



l [m] 



m 2 



400- 



E v [GeV] 



(122) 



m2[eV 2 /c 4 ]' 

The observation of this component together with the neutrino's energy would 
make a determination of the neutrino absolute mass possible. The neutrino's 
energy is measured with uncertainty AE U , which is of the order of 0.1 x E v . 
This uncertainty is 100 [MeV] for the energy 1 [GeV] and is accidentally same 
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Fig. 6: The total probability integrated over the angle per time at a finite 
distance L is given. The constant shows the normal term and the diffraction 
term is written on top of the normal term. The horizontal axis shows the 
distance in [m] and the probability of the normal term is normalized to 0.8. 
Clear excess and decreasing behavior are seen in the distance below 600 [m]. 
The neutrino mass, pion energy, neutrino energy are 1 [eV/c 2 ], 4 [GeV], and 
100 [MeV]. Target is 16 0. 



order as that of the minimum uncertainty h/ Sx derived from the nuclear size 
5x. The total probability for a larger value of energy uncertainty is easily 
computed using Eq. (II 18ft . Figs. 4-7 show the distance dependence of the 
probability. If the mass is around 1 [eV/c 2 ] the excess of the neutrino flux of 
about 20 percent at the distance less than a few hundred meters is found. We 
use mainly m v = 1 [eV/c 2 ] throughout this section. Because the probability 
has a constant term and the T-dependent term, the T-dependent term is 
extracted easily by subtracting the constant term from the total probability. 
The slowly decreasing component decreases with the scale determined by the 
neutrino's mass and the energy. 
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Fig. 7: The probability integrated over the neutrino angle per time at a finite 
distance L is given. The constant shows the normal term and the diffraction 
term is written on top of the normal term. The horizontal axis shows the 
distance in [m] and the probability is normalized to 0.8. Clear excess is seen 
in the distance below 1200 [m]. The neutrino mass, pion energy, neutrino 
energy are 0.1 [eV/c 2 ], 4 [GeV], and 10 [MeV]. Target is 16 0. 



We plot the figure for m v = 0.1 [eV/c 2 ], E v = 10 [MeV] in Fig.0 A 
decreasing behavior is clearly seen. So in order to observe the slowly de- 
creasing behavior for the small neutrino mass less than or about the same 
as 0.1 [eV/c 2 ], the electron neutrino should be used. The decay of the muon 
and others will be studied in a forthcoming paper. 



4.2.3 Energy dependence 

The energy spectrum of the neutrino from the high-energy pion is studied 
next. Since the diffraction term has the origin in the final states that do not 
conserve the kinetic energy, that should show unusual behavior. In Fig. El 
the spectrum for the neutrino mass and pion energy, 1.0 [eV/c 2 ] and 4 [GeV], 
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Fig. 8: The energy dependence of the probability integrated over the angle 
at distance L = 100 [m] is given. The lower curve shows the normal term and 
the diffraction term is added on top of the normal term. The horizontal axis 
shows the neutrino energy in [MeV] and the probability of the normal term 
is normalized to 0.8. The neutrino mass and pion energy are 1.0 [eV/c 2 ] and 
4 [GeV]. Target is 16 0. 

are given. The spectrum of the normal term is flat because the energy in the 
rest system is fixed to one value from the energy-momentum conservation, 
whereas that of the diffraction is not fixed to one value at the rest system and 
is not flat but has a maximum at the energy E v pa E VtJnax /3. The diffraction 
term becomes much larger in much higher energy. 

A unique property of the neutrino diffraction is identified by its energy 
spectrum. The energy spectrum of the normal term from a pion at rest for 
the wave packet size of the momentum width 5 [MeV] is given in Fig. [HI 
The spectrum has a peak at the value derived from the energy-momentum 
conservation, 

m n = E u + p u +p fl = 0, (123) 
of the two body decay. The neutrino energy is uniquely determined to the 
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Fig. 9: The neutrino energy spectrum of the normal and diffraction terms 
in the rest system of the pion are given for the wave packet size of 5 [MeV]. 
The former spectrum becomes wide due to the wave packet effect and the 
latter spectrum becomes wider than the normal component and is in the 
lower energy region. The magnitude of the diffraction term is arbitrary. The 
neutrino mass is 1.0 [eV/c 2 ]. The length is L = 10 [m]. 



value 

2 2 

mt — m„ 

E v = — n * (124) 

The spectrum becomes broad due to the finite wave packet effect. 

The same figure shows another component in the low energy region, which 
corresponds to the diffraction term. In the low energy region of the neutrino, 
the light-cone singularity is not dominant and next to leading terms con- 
tribute. This figure does not include these non-leading terms, hence the 
magnitude of the diffraction term is arbitrary. The length is L = 10 [m]. 

Fig. [TOl shows the energy spectrum of the fraction of the diffraction term 
over the normal term, which are computed with the (V— A) x (V — A) current 
interaction and is represented in a latter section(Sec.7.3) for 56 Fe. Its energy 
is different from Eq. f ll 24j) and spreads over in a wide energy regions. Only 
the leading term is taken into account in this figure. Because the energy 
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Fig. 10: The neutrino energy spectrum of the fraction of the diffraction term 
are given for the energy of the pion 10, 40, 100 [MeV/c] and the length 
L = 10 [m], and for the energy of the pion 60, 100 [MeV/c] and the length 
L = 20 [m]. Target is 56 Fe. The fraction is small in lower energy and is 
larger in higher energy. The diffraction term may be observable in these 
energy regions too. The neutrino mass is 1.0 [eV/c 2 ]. 



is away from the normal term, this component may look like a background 
noise which is uncorrelated with the system. The magnitude becomes small 
in a lower energy region, because the neutrino spreads uniformly in all angle 
there. It is noted that the finite-size correction is not invariant under the 
Lorentz transformation and the magnitude of the diffraction term becomes 
larger in higher energy. Thus the fraction of the electron mode varies with the 
pion's energy. This unusual behavior is a characteristic feature of the 
diffraction component. Our result, in fact, shows that this background 
becomes larger as the pion's energy becomes larger but has the universal 
property. 

4.2.4 Wide distribution of pion momentum 

When a momentum distribution pexpiPn) of initial pions is known, an energy- 
dependent probability is computed using the expression Eq. (11181) . Eq. (11181) 
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is also independent of the position X n and depends upon a pion momentum 
and a neutrino momentum and the time interval T = T„ — TV In experi- 
ments, a position of a pion is not measured and an average over a position 
is made. An average probability agrees with Eq. (11181) . This probability 
varies slowly with the pion's momentum and is regarded constant in the en- 
ergy range of the order of 100 [MeV] . So the experimental observation of the 
diffraction term is quite easy. 

4.3 On the universality of the diffraction term 

The finite-size correction of the probability to detect the neutrino has various 
unique properties. This component is decreasing with the time interval T, 
hence the total probability is not proportional to T in this region. In clas- 
sical particle's decay, the decay process occurs randomly and follow Markov 
process. Hence an average number of decay products is necessary propor- 
tional to T. Now due to the finite-size correction, this property does not hold. 
This is not surprising in the interference region L < l , because the quantum 
mechanical interference effect modifies the probability. 

The finite-size correction is expressed with the universal function g(u u , T), 

2 

where u u = This is determined only with the mass and energy of the 
neutrino and is independent of details of other parameters of the system 
such as the size, shape, and position of the wave packets and others. Hence 
the diffraction component has the genuine property of the wave function 
|muon, neutrino (t)), of Eq. ffl3|) . and is capable of experimental measure- 
ments. 

4.3.1 Violation of conservation laws 

The probability is computed with S[T] and reflects the wave function at a 
finite time t, hence the states of non-conserving kinetic energy contribute to 
the finite-size correction. So conservation laws that are connected with the 
space-time symmetry get modified and various probabilities become different 
from those of T — > oo. The leading finite-size corrections have, nevertheless, 
universal forms that are proportional to g(u, T). 
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4.3.2 Comparisons on the neutrino diffraction with diffraction of 
classical waves through a hole 

1 Inelastic channel. 

The neutrino diffraction is the quantum mechanical phenomenon. The 
neutrino produced in the weak decay of the pion is expressed with the many- 
body wave function composed of the pion, muon, and neutrino. Hence the 
probability to detect the neutrino is computed with this many-body wave 
function. Since this neutrino inside the coherence length is very different 
from the free isolated neutrino, its probability receives the large finite-size 
correction of the universal behavior. Its magnitude, however, depends on 
the wave packet size which is determined with the nucleus that the neutrino 
interacts with. So the finite-size correction is determined by the many-body 
wave function and the out-going wave. We should note that quantum me- 
chanical probability is determined with the overlap of the in-coming state 
with the out-going state and depends on the both states. 

In a classical wave phenomenon, on the other hand, an intensity is deter- 
mined with only the in-coming wave. A magnitude of the in-coming wave is 
directly observed. Hence the finite-size correction and the interference pat- 
tern are determined only by the in-coming wave. Thus interference of the 
quantum mechanical wave is different from that of the classical wave. 

2 Pattern in longitudinal direction 

The neutrino diffraction is a part of the finite-size correction that results 
from the wave natures of the wave function at a finite time t. They are 
generated by the states that are orthogonal to the states at t — > oo hence 
its magnitude is positive semi-definite and depends on the time interval T. 
Hence the neutrino flux has the excess that decreases with the distance in the 
direction to the neutrino momentum and vanishes at the infinite distance. 

The diffraction pattern of light through a hole or the interference pattern 
of light in a double slit experiment are different. The intensity have modu- 
lations in the perpendicular direction to the wave vector. The interference 
term is a product of two waves of different phases and so oscillates. Inte- 
grating the intensity over the whole screen, the oscillating interference terms 
cancel and the total intensity is constant. 

Thus the diffraction pattern of the neutrino is very different from that of 
light. 

3 uj v is Einstein minus de Broglie frequencies 

The pattern of the neutrino diffraction is determined with the angular 
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velocity w v ,diff — ^u,e — ^u,dB- Since oj v ^ e and uj^b are almost the same for 
the neutrino, they are almost cancelled and 0J Vj diff becomes extremely small 
and stable in E v . 

The interference pattern of the light on the screen of the double slit 
experiment, on the other hand, is determined with the angular velocity ou^dB- 
Since co ljdB is large and proportional to the energy, the pattern varies rapidly 
with the position in the screen and with the energy. The diffraction pattern 
of the light passing through the hole varies rapidly. 

4.3.3 Muon in the pion decay 

In experiments of observing the muon in the pion decays, the neutrino is 
not observed. In this situation, the muon's diffraction term has a magnitude 
that is determined by the ratio of the mass and energy, rrvf l /(2E fl ). Since 
the muon mass is larger than the neutrino mass by 10 8 , the value rr? j (2E ft ) 
for the muon is much larger than that of the neutrino by 10 16 . For the 
muon of energy 1 [GeV], the length is of the order of Iq = 10 -14 [m]. This 
value is a microscopic size and p(T,u; M ) vanishes at a macroscopic length. 
Hence the probability of detecting the muon at the macroscopic distance 
becomes constant. The muon from the pion decay shows neither the finite- 
size correction nor the diffraction effect. This probability agrees with the 
production probability. The muon and neutrino behave differently at the 
finite distance. 

If the muon is observed under a condition that the neutrino is detected at 
the finite T, S[T] is applied and the probability to detect the muon has the 
contribution from the neutrino diffraction. The diffraction component gives 
a wide energy spectrum for the muon since that comes from the tail of the 
distribution function. Fig. [TT1 shows a probability integrated over the neutrino 
energy in this condition that both the muon and neutrino are detected, which 
is obtained from Eq. (|4"H|). In this figure, we use units c = l,h— 1 and express 
the energy and time with the neutrino mass m u . Energy of the pion is 250 
m v and the muon has the energy 210 m v and has an angle with the pion of 
cos 9 = 0.95 — 1. The cosine of the angle between p n — and p v is in the 
horizontal axis. T is 30 m" 1 . The neutrino mass of an unphysical magnitude 
of the order of MeV and the value of T are chosen in such manner that the 
numerical computation of diffraction component is easily made. Qualitative 
features of Fig. [11] are that there exist a large peak at cos 9 ~ 1 and small 
peaks at the tail region. The former is the peak from the root of u = at 
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Fig. 11: The angle dependence of the diffraction and normal components of 
probability is given. The large peak shows the normal component and the 
small peaks at the tail of the previous peak show the diffraction component. 
The horizontal axis shows the cosine of the angle between p v and p n — p^ and 
the vertical axis shows the probability. The pion energy, muon energy, and 
time interval are 250 m u , 210 m v , and 30 m" 1 . 

bp ~ and 5E = and the latter are the peaks from the roots of u = of 
5p 7^ and SE ^ 0. The latter peaks, which do not exist in the probability of 
detecting only the muon, show the feature of the diffraction component of the 
probability when the neutrino is detected. Thus the diffraction component 
is observed in the muon also when the neutrino is detected simultaneously. 
Experimental verification of the diffraction term of this situation using the 
muon may be made in future. 

As was shown in Section 3, the production rate is common to the muon 
and neutrino, since they are produced in the same decay process. However the 
rates measured by apparatus depend on the condition of the measurement. 
If one particle is measured and other is un-measured, its rate for the neutrino 
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receives the large finite-size correction, but the rate for the muon receives no 
correction. Hence they become different each other. If both particles are 
measured simultaneously, the rates for both are the same. 

Thus the neutrino is in the non-asymptotic region of the finite correlation 
in wide area, and the transition amplitude for a neutrino that is observed at 
a finite distance by a nucleus becomes different from that of the infinite dis- 
tance. The neutrino wave is a superposition of those waves that are produced 
at different positions and the probability gets an additional contribution and 
the probability is modified by the diffraction term. The overlap between 
the neutrino wave that is detected with a nucleus in a detector and those 
that are produced from a pion decay shows the neutrino diffraction of unique 
properties. So the neutrino flux measured with its collisions with a nucleus 
in targets are different from that defined from the norm of wave function. 



5 Implications 

In this section, various physical quantities of neutrino processes which are 
modified by the neutrino diffraction are studied. Particularly neutrino nu- 
cleon total cross sections, quasi-elastic cross sections, electron-neutrino pro- 
duction anomaly, a proton target enhancement, and an anomaly in atmo- 
spheric neutrino are such processes that have significant contributions from 
the neutrino diffraction. 



5.1 Total cross sections of z^-N scattering 

Neutrino collisions with hadrons in high energy regions are understood well 
with the quark-parton model. A total cross section of a high energy neutrino 
is proportional to the energy and is written in the form 

C = ¥Z^(Q + 1/3Q), (125) 

7T 

using integrals of quark-parton distribution functions q(x) and q(x) and Q = 
dxxq(x), Q = Jq 1 dxxq(x). The cross section is proportional to the neutrino 

energy and a current value is a u /E = 0.67 x lCr 38 [cm 2 /GeV]. 

Now the rate of process of the neutrino produced by a decay of a pion 

and interacting with a nucleus has a finite-size correction. It modifies the 

probability of the neutrino collision and the cross section. We estimate its 
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effect hereafter. Including the diffraction term, the effective neutrino flux 
becomes the sum of the normal and diffraction terms 

/ — fnormal(l + r diff), (126) 

where rjiff is the rate of the diffraction component over the normal compo- 

2 

nent, and is a function of the combination (^-L), 

2 

TYt 

rdlff = d °~ 9i 2^k L) > (127) 

where L is a length of the decay volume and the coefficient do is determined 
from geometries of experiments. 

When a detector is located at the end of the decay volume, the correction 
factor Eq. (11 2 7p is used. In actual case, the detector is located in a distant 
region from the decay volume. There are material or soil between them and 
pions are stopped in beam dump. The neutrino is produced in the decay 
region and propagates freely afterward. Since the wave packets of one a v 
form the complete set [10J, the wave packet of the size at the decay volume 
is the a v determined with the detector. The neutrino flux at the end of the 
decay volume is computed with the diffraction term of the decay volume 's 
length L and the wave packet size of the detector. Wave packets of this a v 
propagate freely from the end of decay volume to the detector. The final 
value of neutrino flux at the detector is found combining both effects. When 
neutrino changes flavor in this period, the final probability for each flavor is 
written with a usual formula of flavor oscillation. 

The true neutrino events in experiment is converted to the cross section, 
a exp (E), that includes the diffraction component and is connected with the 
cross section computed with only the normal component a the (E) by the rate 

a the (E) = a exp (E) . 

1 + ^// 

Conversely the experimental cross section is written as 

a exp (E)/E = (1 + r di ff)(a the (E)/E). (128) 

a the (E)/E is constant from Eq. (11251) so the E-dependence of a(E) exp / E is 
due to E-dependence of r<&//, Eq. (I127p . 



48 



Global average 

MINOS (DATA) = A - 
NOMAD(DATA) -{3-' 




Neutrino energy [GeV] 

(a) 



MINOS(diffraction) Q 
NOMAD(Diffraction) ^ 
Global average 



O 

X g 

x rv 
X Q 

X, 




Neutrino energy [GeV] 
(b) 



Fig. 12: Neutrino- Nucleon total cross section of MINOS and NOMAD (a) and 
total cross sections of the sums of normal and diffraction terms in geometries 
of MINOS and NOMAD (b) are given. The horizontal axis shows the neutrino 
energies in [GeV] and the vertical axis shows the ratio of the cross section 
over the energy. 



The correction r^// depends on the geometry of experiments and the 
material of the detector. We compute r di ff using the experimental conditions 
of MINOS [30] and NOMAD [3T] and the total cross sections. The geometry 
of MINOS and NOMAD are the following. The lengths between the pion 
source and the neutrino detector, L det _ so , and those of the decay region, 

^ decay— t eg i are. 

NOMAD : L det „ so = 835 [m], L decay _ reg = 290 [m], (129) 
MINOS : L det _ so = 1040 [m], L decay _ reg = 675 [m]. (130) 

Also pion beam spreading was included from angle of initial pion; to 
10 [mrad] for NOMAD and to 15 [mrad] for MINOS. 

The wave packet size is estimated with the size of target nucleus. From 
the size of the nucleus of the mass number A, we have a v = As /m^. For 
various material the value are 

o v = 5.2/m^; 12 C nucleus. 

o v = 14.6/to^; 56 Fe nucleus. (131) 

Including the geometries, beam spreadings, and wave packet sizes, we 
computed the total cross sections and compared with the experiments in 
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Fig. [T21 These cross sections computed theoretically slowly decrease with 
the energy in the geometry dependent manner and agree with the experi- 
ments. Since the experimental parameters such as the neutrino energy and 
others are different in two experiments, the agreements of the theory with 
the experiments are highly non-trivial. So the large cross sections at low 
energy regions may be attributed to the diffraction component. 

We have compared only NOMAD and MINOS here. Many experiments 
are listed in particle data [32] and most of them have similar energy depen- 
dences and agree qualitatively with the presence of the diffraction compo- 
nents. It is important to notice that the magnitude of diffraction component 
is sensitive to geometry. Furthermore if a kinematical constraint Eq. flllOp on 
the angle between p n and p u was required, only the events of the normal term 
was selected. Then the cross section should agree with that of the normal 
term. 

5.2 Quasi-elastic cross sections 

Quasi-elastic or one pion production processes are understood relatively well 
theoretically. The diffraction modifies the total events of these processes also. 
The cross sections for 

v + n^ //" +p(+vr°), (132) 

V+P -+ H~ +P + 7T+, (133) 

u + p^ fi + + n(+7r°), (134) 

and the neutral current process 

u + N ->■ v + N(+7r°), (135) 

are known well using CVC, PCAC, and vector dominance and are studied 
recently by MiniBooNE [33] • The parameter is the axial vector meson Ma 
and higher mass contributions. So these cross section are used to study the 
diffraction terms. 

5.3 Electron neutrino anomaly 

In pion decays, a branching ratio of an electron mode is smaller than that 
of a muon mode by factor 10~ 4 due to the helicity suppression of the decay 
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of a pseudo-scalar particle caused by the charged current interaction. This 
behavior of the total rates has been confirmed by the observations of charged 
leptons. 

Now the probability to detect a neutrino inside the coherence length, 
where the neutrino retains the wave natures, is affected by the finite-size 
correction. Because this correction comes from the states that have different 
kinetic-energy from the initial value, the neutrino in this region does not 
follow the conservation law satisfied in the asymptotic region t — > oo. The 
rate that electron neutrino is detected is not suppressed. The ratio of the 
probability to detect the electron neutrino over that of the muon neutrino 
becomes substantially larger in near-detector regions. 

To compute the transition probability and the spectra of electron and 
muon neutrinos, we start from the (V — A) x (V — A) interaction Lagrangian 
(Hamiltonian). The result of the probability is almost the same as that 
of Eq. (jl]) in the muon mode but is different in the electron mode since 
the diffraction component does not satisfy the rigorous conservation of the 
kinetic-energy and momentum. In I, it was found that the initial pion is 
described by a wave packet of a large size. Hence the initial pion of the plane 
wave is studied here. The amplitude T is written with the hadronic V — A 
current and Dirac spinors in the form 



where N = ig(a u /7i) Ti {mim u / EiE u ) 2 , and the time t is integrated in the 
region T^ < t < T u . The transition probability to this final state is written, 
after the spin summations are made, with the correlation function and the 
neutrino wave function in the form 



where C = g 2 (Att jo v ) 2 V 1 , V is a normalization volume for the initial pion, 



T 




(136) 



4 




3 




(2k) 




(2(Ptt • Pu){Pn ■ Pi) - ml 



(pi-Pu))e-^- pi > 



(138) 
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The probability of detecting a neutrino of p u at X v and a lepton I of 
arbitrary momentum is expressed as the sum of the normal term Go and the 
diffraction term g(T,u u ), 

P=N2 f g ? fVfvK g ; 2fVfV) BCT |M J + Go] , (139) 

where ./V2 = 8Tg 2 ov and L = cT is the length of decay region. In Go the 
energy and momentum are conserved approximately well and 

Pi~Pn- Pu, (140) 

is satisfied. Hence from a square of the both hand sides, the mass shell 
condition 

rrvf « m 2 - 2p n ■ p v , (141) 

is obtained. Thus the normal terms are proportional to the square of lepton 
masses and the electron mode is suppressed pEJ El El [9] . In g(T,u v ), on the 
other hand, momenta satisfy 

Pl^Pn- Pu, (142) 

and the diffraction terms are not proportional to the square of lepton masses 
and the electron mode is not suppressed. 

The total probability of detecting a neutrino or a charged lepton in the 
pion decay at macroscopic distance is written in the form, 

P = Pnormal + Pdiff- (143) 

In Eq. (I143p . P noTma i is the normal term that is obtained from the decay 
probability Go in Eq. f ll39[) and -P<t/7 is t ne diffraction term that is determined 
from g in Eq. ( 1 1 3 9 [) . The former probability agrees to that obtained using 
the plane waves and the latter one has not been included before and its effect 
is estimated here. The diffraction term at T is described with its mass and 
energy in the universal form 

li ff = CTg{T, Ul ). (144) 

The frequency ui are small for neutrinos and large in charged leptons. g(T, uj{) 
is positive definite and decreases slowly with a distance L and vanishes at 
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infinite distance. Hence at L = oo, the probability agrees with the normal 
component, 

P Pnormal- (^-^v 

The length scale 2cE[/mf is macroscopic size in neutrinos but is 10 -10 [m] or 
less for the electron and muon. The magnitude of g(T, ooi) at the macroscopic 
distance is given in Fig. 2 of I. At L = 100 [m], E = 1 [GeV] for the mass 
1 [eV/c 2 ] (v), 0.5 [MeV/c 2 ] (e), and 100 [MeV/c 2 ] (ji), the values are, 

g(T,u v ) w 3, 

S(i> e ) «o, 

j(T,w^0. (146) 

In this region, they satisfy 

2 

S(T,o;,) « -iSClNa;,,), (147) 

hence the diffraction component at a macroscopic distance is finite in the 
neutrino and vanishes in others. It is striking that the probability to detect 
the neutrino has an additional term and is not equivalent to that of the 
charged lepton even though they are produced in the same decay process. 

The diffraction term is generated by the tiny neutrino mass and the light- 
cone singularity. Hence the pattern is determined by (E(p u ) — c\p u \)t and 
becomes long-range. Because Pdiff is the finite-size correction caused by the 
neutrino interference, it has different properties from P nor mai in flavor and 
momentum dependences. The neutrino diffraction furthermore is sensitive 
to the absolute neutrino mass. We study implications of the diffraction term 
to the electron mode in the pion decays here. In Fig.[13l experiments of 
LSND [31] and the two neutrino experiment ( TWN) [35] are compared with 
the diffraction components and the flavor oscillations. Theoretical values 
are obtained including geometries of the experiments. Since those of LSND 
and TWN are different, the theoretical value for the LSND are smaller than 
that for TWN. The experimental values plotted with crosses agree with the 
theoretical values. The values from the flavor oscillations expected from the 
current parameters are also shown. The mass-squared differences and mixing 
angles from the recent ground experiments lead negligible values for both 
experiments. A sterile neutrino of the mass around 1 [eV/c 2 ] is necessary 
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Fig. 13: Experiments of LSND and TWN are compared with the theoreti- 
cal values of the diffraction terms. TWN(EXP) and LSND(EXP) show the 
experimental values and TWN (Diffraction) is computed with the parame- 
ters m v = 0.2[eV/c 2 ], E v = 250[MeV], P n = 2[GeV/c], LSND (Diffraction) 
is computed with m u = 0.2 [eV/c 2 ], E v = 60[MeV], P n = 300[MeV/c]. 
Flavor oscillation oscillation(T2K) shows the values for sin 2 # 13 = 0.11, 
5m 2 23 = 2.4 x 10~ 3 [eV 2 /c 4 ], E v = 60 [MeV], and LSND(sterile) shows with 
sin 2 fl = 0.004, 5m 2 = 1.2 [eV 2 /c 4 ], E v = 60 [MeV]. 

to fit the data of LSND with the flavor oscillation. The agreements of the 
values from the neutrino diffraction in LSND and TWN suggest that it is 
unnecessary to introduce additional parameters. 

In Fig. HI], the maximum possible fraction of the electron neutrino in a 
geometry of T2K experiment is shown. The spreading of the pion beam is 
ignored in this Figure. Since the diffraction is sensitive to the pion beam 
spreading, the real value may becomes smaller than this figure. 
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Fig. 14: Fraction of the electron neutrino of the mass 0.2 [eV/c 2 ], 0.5 [eV/c 2 ] 
and l[eV/c 2 ] at L=110 [m], distance=170 [m] of T2K geometry and P n = 2 
[GeV/c]. 

5.4 Proton target anomaly 

Magnitude of diffraction component depends upon the size of the nucleus 
which neutrino interacts with and is expressed by the wave packet size o v . 
It becomes larger with the larger target. It is known and used in the text 

2 

that nuclear size is proportional to As and the large A nuclear gives a large 
diffraction component, generally. Proton has a smallest intrinsic size. How- 
ever a proton is expressed by a wave function of its position in matter. So the 
wave packet size is determined by a size of this wave function. Since a proton 
is the lightest nucleus, it has the largest size. We estimate this size using 
center of mass gravity effect between proton and electron. For the proton's 
mass m p and the electron mass m e , an electron's coordinate x e i ec t ron and the 
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proton coordinate expressed as, 

Xelectron = X H — f» X + (l + 7—^—7)^, 

m e + m p 2000 

x p = X — r^X — r. (148) 

p m e + m p 2000 v ; 

If the wave function of the atom is 

tt(flMf), (149) 

and the function of the relative coordinate, tp(r), is extended by an amount 
Ratom which is about 10~ 10 [m] then the proton is extended with a radius 

R P = 2m R atom ~ 5 x 10~ 14 [m]. (150) 

This value is much shorter than the atomic scale and is larger than one 
nucleon's size 1 [fm] = 10~ 15 [m] by factor 50. Thus proton in solid is extended 
to the size of the atomic wave function, which is larger than the nuclear 
size of O. Hence proton gives the important role in the neutrino diffraction. 
Its size may be in the range 

lproton(U) = 5 x 10- 14 - 10- 13 [m]. (151) 

An enhancement of diffraction contribution due to the proton is expected 

in 

u + p -> /i + + X°, X° = n, pir', mr°, others. (152) 



5.5 Atmospheric neutrino 

The neutrino flavor oscillation was found first with an atmospheric neutrino. 
Neutrinos are produced from decays of charged pions and muons in secondary 
cosmic rays. Since the matter density is low in atmosphere, these charged par- 
ticles travel freely long distance. Thus neutrinos produced in decays of pions 
or muons show the diffraction phenomenon and the diffraction components 
are added to the neutrino fluxes. These neutrino events may be observed in 
detectors set in the ground, such as Super-KamiokaNDE(SK) if the absolute 
mass is a reasonable value. The minimum mass allowed from the mass- 
squared difference is about the value, V 5m? 10 -2 [eV/c 2 ]. Then the length 
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that the diffraction component is observed becomes L = 20 [km] for 

E v — 1 [GeV], which is longer than the height of troposphere. Hence the 
diffraction component could be observed with the angle-dependent excess 
of the electron and muon neutrino fluxes. Since the diffraction components 
from pion decays are common to both neutrinos, their ratio is not sensitive 
to the diffraction. Instead of this ratio, a ratio of the neutrino flux to the 
flux of charged leptons is good to see the signal of the neutrino diffraction. 

6 Behaviors of the probability suggest vio- 
lations at first sight but do not so in fact. 

6.1 Unitarity 

Probability of detecting the neutrino per time -P(L) decreases with the dis- 
tance L. This behavior of the probability appears to suggest that the prob- 
ability is not preserved and is inconsistent with the unitarity. However this 
behavior is derived from the S[T] that satisfies S , '[T] J S'[T] = 1 and is consis- 
tent with the unitarity. The probability at L is determined with S-matrix 
S[T], L = cT and has two components P = p("™0 + P di ff(h). Both 
terms are positive semi-definite and the latter is decreasing with L because 
the constraint to the final state from the energy conservation becomes more 
stringent with increasing L. This decreasing behavior is a natural consequence 
of the unusual properties of the finite-size correction and is consistent with 
the unitarity S , t[T]S'[T] = 1. The unitarity leads that the life time of the 
pion becomes larger if the neutrino is detected at a finite T. 

6.2 Lepton number non- conservation 

The probability of the pion decay process in the situation where the neutrino 
is detected has the large finite-size correction. The charged lepton shows the 
same behavior in the same situation. This will be confirmed if the charged 
lepton is measured simultaneously with the neutrino in experiments. This 
has not been done and is consistent with the lepton number conservation. 

The probability of the pion decay process in another situation where the 
neutrino is un-detected has no finite-size correction. The charged lepton in 
this process does not show the finite-size correction and its decay probability 
is computed with the standard calculation of using the plane waves. This 
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situation has been studied well experimentally and agrees with the theoretical 
calculations obtained with £[00]. 

Now the boundary conditions of the above two cases are different. One 
boundary condition leads uniquely one consequence and the different bound- 
ary conditions may lead the different probabilities. Our results show that 
the different boundary condition on the neutrino leads the different result 
on the decay probability. Thus the probability to detect the neutrino in the 
first case is different from that of the charged lepton in the second case. It is 
meaningless to compare the probability for neutrino in the first case with that 
for the charged lepton in the second case, because they follow the different 
boundary conditions. 

Since the probability of detecting the neutrino at a finite distance deviates 
from that at the infinity and the probability of producing the neutrino is de- 
fined with the value at the infinite distance, they are different. The fact that 
the probability of detecting the neutrino is different from that of the charged 
lepton does not mean the violation of the lepton number conservation, but 
means that the probabilities depend on the boundary condition. The lepton 
number is conserved. Thus the different behavior of the probability from 
that of the production is similar to that of the retarded electric potential of 
a moving charged body. 

6.3 Dependence on wave packet size 

It is known that the total probability at T = 00 does not depend on the 
wave packet size [21] . The result of the present paper Eq. (If 18|) in fact shows 
that the first term in the right-hand side is independent of the wave packet 
size. Now the second term in Eq. (If f 8p . which is the finite-size correction, 
is proportional to a v . Since S[T] is determined with the wave packets, its 
boundary conditions are determined with <j v and the finite-size correction 
depends on o v . That increases with o v and diverges at <j v = 00. The 
diverging correction at a u = 00 is consistent, in fact, with the fact that 
the total cross section diverges for the plane waves [21] . This occurs because 
the denominator of the neutrino propagator vanishes. Nevertheless the finite- 
size correction has the universal properties, which were computed with S[T] 
defined by the wave packets. 
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7 Unusual features 



7.1 Energy non-conservation and violation of symme- 
tries 

The S-matrix at the finite-time interval S[T] does not commute with the free 
Hamiltonian Hq, and satisfies Eq. (I22j) . In this region, one pion and decay 
products co-exist in a coherent manner that is determined by Hi. Hence 
Hi has a finite expectation value and the kinetic energy defined by the free 
part H is different from that of H. H is not conserved in S[T], despite the 
fact that the total energy defined by H is conserved. The contribution to the 
transition probability from these states was computed in the text analytically. 
They exhibit the neutrino diffraction and give the finite-size correction to the 
neutrino flux. 

The finite-size correction is not invariant under Lorenz transformation 
either and the magnitude of the diffraction term in the pion rest system is 
smaller than that in the high energy pion. 

7.2 Helicity suppression 

Decay rate of the pion to the electron mode is suppressed over that to the 
muon mode by the helicity suppression. The helicity suppression hold in a 
decay of a pseudo-scalar particle to a neutrino and lepton caused by V — A 
weak interaction. Conservations of the kinetic energy and the angular mo- 
mentum enforce vanishing of the amplitude for the massless lepton. Now the 
kinetic energy is not conserved in the non-asymptotic region and the helicity 
suppression is not effective. The probabilities to detect the neutrino in this 
region are not suppressed. Thus the normal terms hold these properties and 
the electron mode is suppressed, whereas the finite-size correction does not 
hold these properties. The electron mode is not suppressed in the diffraction 
component and has the sizable magnitude in the spatial region where the 
normal term is negligible. 

Thus when the neutrino is observed in near-detector region, the electron 
neutrino is substantially enhanced. 
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7.3 Large finite-size correction 

The finite-size correction of the probability to detect the neutrino is finite in 
the spatial region L < / of the distance L between the pion and the neutrino, 
where the interaction energy is finite and the energy defined by the free 
Hamiltonian H deviates from the conserved total-energy. The states of these 
energies of H Q that are different from that of the initial state interact with 
S[T] in a universal manner that is determined by Hi and gave the finite-size 
correction. Particle spectrum at ultra-violet region is universal in relativistic 
invariant systems and gives the light-cone singularity to correlation functions. 
The light-cone singularity is real and extended to a large area and gives the 
universal finite-size correction to light particles. It might be remarkable that 
the states of the ultraviolet region give the observable effect to the probability 
of the tree diagram. 

Since the probability to detect the neutrino has the large finite-size cor- 
rection, the neutrino detection affects strongly the pion decay probability. 
The decay rate of the pion in the situation where the neutrino is detected 
is different from that where the neutrino is un-detected, because they are 
described by the amplitude of the different boundary conditions. 

The life time of the pion in which the produced neutrino is observed 
becomes shorter than the normal value. This phenomenon that the life time is 
modified by its interaction with matters is known in the literature as quantum 
Zeno effect. Neutrinos actually interact extremely weakly with matters and 
a majority of neutrinos are passing freely without any interaction and are 
not affected by this effect. Consequently the majority of the pions are not 
affected by the finite-size effect and its life time is not modified and has the 
normal life time. Although the detected neutrino receives the large finite-size 
correction, its effect is negligibly small for observables of the pions. 

7.4 Interference in parallel direction to the momentum 

The diffraction term depends on the time interval T and is positive defi- 
nite. Its pattern changes with the distance L extremely slowly. Thus the 
interference pattern is along the parallel direction to the momentum. 
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7.5 Dependence on the apparatus 

The neutrino diffraction appears in the quantum mechanical transition prob- 
ability and depends on both the initial and final states. The wave function 
of the final state is determined by the apparatus and so are the scattering 
amplitude and probability. It is in fact quite reasonable that the interfer- 
ence pattern depends on the apparatus in quantum mechanics. In classical 
physics, on the other hand, any physical variables are observable and in- 
terference patterns do not depend on the apparatus. Hence the fact that 
interference pattern of the probability of the present work depends upon 
the apparatus and is quite different from that of the classical physics is not 
surprising. 

8 Summary and implications 

We found that the neutrino in the pion decay has a large non-asymptotic 
region of having wave natures and the probability to detect the neutrino in 
this region reveals unusual property. The neutrino is unique and shows the 
large finite-size correction of the diffraction form. The unusual diffraction 
pattern is caused by the relativistic invariance and the tiny neutrino mass. 
Its origin, mechanism, characteristic features, and implications are presented. 

The pion, neutrino, and charged lepton are described by a many-body 
wave function that follows the Schrodinger equation. The wave function 
at a finite t has various unusual properties that are different from those at 
t — > oo. Especially the total kinetic energy is different from that of the 
initial state and takes a wide range of values. Hence the wave functions of 
these particles vary in space and time. The physical quantities of these wave 
functions are computed by the S-matrix at a finite-time interval that satisfies 
the boundary condition at a finite time. The position-dependent probability 
to detect neutrino, which is proportional to the flux, was computed with the 
S[T] expressed by wave packets and the large finite-size correction of universal 
property that is insensitive to the pion's initial conditions was found. The 
pion's mean free path estimated in Appendix of I shows that this is long 
enough for the new term to be observed in the experiments. The diffraction 
term reflects the finite-size correction, and depends on the distance, energy, 
angle, and the absolute neutrino mass in the universal manner. 

The wave function at a finite t has the continuous kinetic energy that 
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leads the space-time dependent wave nature, which is probed with the wave 
packet localized in space. Eq. (148 j) is derived in this manner and shows a 
reason why the probability has the large finite-size correction. The angular 
velocity u is given by the energy of the moving frame u = 5E — v v ■ 8p, 
consequently uj = has the root of SE ^ and 5p ^ in addition to that 
of SE = and 5p = 0. The former root exists only in S[T] and gives the 
finite-size correction. The probability around the former root receives also 
the finite-size correction, as is shown also in Appendix B. Since the kinetic 
energy and momentum are not conserved in this root, the finite-size correc- 
tion has various unusual properties. Thus the states which violate the energy 
and momentum conservation give the correction. To compute the dominant 
contribution in the long-distance region rigorously, it is convenient to inter- 
change the order of the integrations and introduce the correlation function 
^^(Sx). The states of non- conserving the energy and momentum form the 
light-cone singularity of A nyfJ- (5x) and give the most important contribution 
to the finite-size correction of the probability to detect the neutrino. The 
regular terms of A n>lx (Sx) give the constant value of the probability. 

The light-cone singularity is formed from the superposition of waves of 
the infinite momentum, which have always the light velocity and vanishing 
complex phase. So when it is combined with the slow phase C of the neutrino 
wave function, Eq. (H0|) . this slow phase remains to physical observable. 

The slow phase of the neutrino wave packet is the outcome of the can- 
cellation between the time and space parts of extremely light neutrino. The 
phase is determined by the difference of space-time coordinates (5t, 5x) and 
the central values of the energy and momentum, as E(p)5t — p-Sx, where 
the energy E(p) is given by y/p 2 + m 2 . When the difference of positions 
Sx is moving with the light velocity in the parallel direction to the mo- 
mentum p, v c = c(p/\p\), as 5x = v c 5t, then the total phase becomes 
E v St — p u -5x = (m 2 ,/2E 1/ )8t. The angular velocity becomes the small value 
m 1 ,/ (2E V ) and makes the interference phenomenon long-range. The new term 
in the probability decreases slowly with the distance in the universal man- 
ner determined by the mass and energy of the neutrino as m 2 J /(2E u ). This 
form is independent from the details of wave packet shape and parameters 
as far as the reality of the neutrino wave function w{x) is satisfied, which is 
ensured from the invariance under the time inversion. The relative magni- 
tude of this component is not universal and depends upon the size of wave 
packet. Based on the estimation of the size, we found that the magnitude of 
the new universal term is sizable for the measurement. Since the slope of the 
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diffraction component is determined by the mass and energy of the neutrino, 
the absolute value of the neutrino mass would be found from the neutrino 
interference experiments. 

The diffraction components were compared with several previous neu- 
trino experiments in Section 7. First, the slight energy dependence of the 
total neutrino-nucleon cross sections at high energy regions, which is hard 
to understand in the standard theory, was shown to agree with the excess 
of the effective neutrino flux due to the diffraction. The excesses of neu- 
trino events will be observed in other reactions as well at macroscopic short 
distance regions. Theoretical calculations at distances of the order of a few 
hundred meters were computed and shown in Figures of Section 6. From 
these figures, the excesses are not large but are sizable magnitudes. Hence 
these excesses shall be observed in these distances. Actually fluxes measured 
in the near detectors of the long-baseline experiments of K2K (36] and Mini- 
BooNE [37] may show excesses of about 10 — 20 percent of the Monte Carlo 
estimations. Monte Carlo estimations of the fluxes are obtained using naive 
decay probabilities and do not have the coherence effects we presented in the 
present work. So the excess of these experiments may be related with the 
excesses due to interferences. The excess is not clear in MINOS [38]. With 
more statistics, qualitative analysis might become possible to test the new 
universal term on the neutrino flux at the finite distance. If the mass is in the 
range from 0.1 [eV/c 2 ] to 2 [eV/c 2 ], the near detectors at T2K, MiniBooNE, 
MINOS and other experiments might be able to measure these signatures. 

Second, the suppression of the electron mode in pion decay is modified in 
the probability to detect neutrinos. Since the energy-momentum conserva- 
tion does not hold in the S[T], the helicity suppression mechanism does not 
work in the diffraction component. So the electron mode is enhanced drasti- 
cally. The theoretical value of fraction that includes the finite-size correction 
of the electron mode was compared with LSND and TWN, and agreements 
were obtained. Further confirmation of the diffraction component by observ- 
ing the electron neutrino in pion decay will be made using modern version of 
LSND or similar experiments. T2K near detector is a possible place for that. 
Third, anomalies in proton target, and atmospheric neutrino would supply 
also specific signature of the neutrino diffraction. The neutrino diffraction 
is sensitive to the absolute neutrino mass but is not so to other parameters 
such as pion and neutrino energies. Hence the observations of the neutrino 
diffraction is easy. The absolute neutrino mass could be found with these 
experiments. 



63 



At the end, we summarize the reasons why the interference term of the 
long-distance behavior emerges in the pion decay and is computed with the 
wave packet representation. The connection of the long-distance interference 
phenomenon of the neutrino with the Heisenberg's uncertainty relation is 
also addressed. 

Relativistic invariance forces the particle's momentum unlimited and makes 
the correlation function A n>fl (St, Sx) have a singularity near A = 0, which is 
extended to large distance \Sx\ — > oo. This is one of the features of rela- 
tivistic quantum fields in Minkowski space-time and is one reason why the 
long-range correlation emerged. For a non-relativistic system, on the other 
hand, the same calculation for stationary states is made by, 

J dk{x!\k){k\x 2 ) = S(x! - x 2 ), (153) 

and the only one point Sx = satisfies the condition. Long-range correlation 
is not generated. The rotational invariant three-dimensional space is com- 
pact but the Lorentz invariant four-dimensional space is non-compact. So 
it is quite natural for the non-relativistic system not to have the long-range 
correlation that the relativistic system has. The light-cone singularity is the 
peculiar property of the relativistic system. 

Heisenberg uncertainty relation is slightly modified in the wave along the 
light cone. The neutrino wave function behaves at large distance along the 
light-cone region in the form 

i(E v t-p„-x) p 1 ^ 1 

Mt,x) = f = /— — , (154) 

Ob Cv 

where / has no dependence on the distance \x\. Consequently the uncertainty 
relation between the energy width SE and the time interval St becomes 

St6^ = K8Ex±(^y*h. (155) 

The ratio [m v j ' E v ) 2 is of the order of 10~ 18 and St becomes macroscopic even 
if the energy width SE is microscopic of the order of 100 [MeV]. For instance 
if the pion Compton wave length, A^, is used for the microscopic length, then 
cSt becomes 

10 18 x A x « 10 3 [m], (156) 
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which is about the distance between the pion source and the near detector 
in fact. So interference effect of the present paper appears in this distance 
and is observable using the apparatus of much smaller size. 

In the time interval T < Iq/c, the finite-size correction is not negligible. 
The diffraction term has the finite value and contributes to the probability. 
Hence, the probability to detect the neutrino deviates from the production 
probability. In another region l /c < T, the neutrino behaves like a free 
isolated particle and is in asymptotic region. The diffraction term vanishes 
and the probability agrees with the normal term. The probability to detect 
the neutrino is computable with the ordinary S-matrix and agrees with the 
production probability. Wave packet formalism is applicable to both of the 
asymptotic and non-asymptotic regions. 

The characteristic small phase of the neutrino wave function of the an- 

2 

gular velocity ^jf- along the light cone causes the diffraction phenomenon 
at the anomalously large area. There would be similar phenomena in other 
light particles or others where the scattering matrix of the finite-time inter- 
val S[T] are important. Unique properties of S[T] may give new insights to 
those phenomena that are hard to understand with S'foo]. 

In this paper we studied the amplitude and probability in the lowest order 
in Gf and ignored higher-order effects such as pion life time, pion mean free 
path, and effects of electroweak gauge theory in studying the quantum effects. 
They do not modify the effect of the tree diagrams studied in the present 
work. We will study these problems and other large scale physical phenomena 
of low energy neutrinos in subsequent papers. 
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Appendix A Long-range correlation for gen- 
eral wave packets 

Non-Gaussian wave packets were studied in the general manner and the uni- 
versal behavior of the phase was obtained in the text. In this appendix, 
the explicit forms of the wave packets are studied as concrete examples. It 
is re-confirmed that the long-range component of the probability at around 
t = 1nE v jvr? v becomes the universal form, 
type 1 



68 



One way to express the non-Gaussian wave packet is to multiply Hermi- 
tian polynomials and to write the amplitude in the form 



(2vT 



where H n is assumed to be real in order for the wave packets to preserve the 
time reversal symmetry and an even function of k v — p u in order for the wave 
packets to preserve parity, as was shown in Appendix of I 
For the simplest case 

H n = a v (k v -p v ) 2 } (158) 
the wave packet in the coordinate representation is 



3 

"J \ n,, 



(159) 

and is substituted into the integral Eq. ( 1681) . 



J SW = NlMi^ E -^[- 1 ^ + JL(1 - v u )\5tf + 0(1 - v v f{5t)% 

(160) 



' u 



Thus the phase factor has the same universal form as the Gaussian wave 
packet and the correction is determined by the negligible small parameter 

^(i-^)W=(^;) 2 (^) 2 - 

We have proved that the correlation function of the non-Gaussian wave 
packet has the same slow phase and long-range term as the Gaussian wave 
packet and the small correction becomes negligible for the simplest case 
Eq. 11581 Hence for any polynomials H n that are invariant under the time 
and space inversions, the correlation function has the same long-range term 
and small negligible corrections. 

type 2 

For another non-Gaussian wave packet 

-^3 [ d Z ve -«&ME&)(t-T v )-UZ-Xu)) (161) 

(2tt)§ J 
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we have the same result, 
type 3 

In the type 1 and 2 the time reversal and parity symmetries are assumed 
for the wave packet shape. If these symmetries are not required, the function 
H n or a has an imaginary part. In this case, the correlation function has a 
correction term in the order of (1 — v)(t\ — t 2 ) and this term is expressed 

(1 - v v )St = ^r^St, (162) 

hence the correction term vanishes at the high energy. With a suitable param- 
eter, the universal form of the slowly decreasing component of the probability 
of the present work may become observable even in arbitrary system. The 
Lorentz invariant form of the energy dependent phase of the wave packet 
and the light-cone singularity of the pion and muon decay vertex give this 
universal behavior. 



Appendix B Origin of the large finite-size cor- 
rection 

We study the reason why the finite-size correction are large for the light par- 
ticles from the amplitude. Plane waves are assumed for the initial pion and 
un-detected particle and the wave packet of a v is assumed for the neutrino 
in the final state. Integrating the coordinate x in Eq. (T28]) . for = oo, we 
have the amplitude, 

T 

iuit 

) 



r = Ce < * «(p^(l-75)«(p,)(0|^_ j4 (0)|7r)e-^ / die 

Jo 

uj = 5E -v v - Sp, 5E = E(p n ) - E(p fJ ) - E{p v ), 5$=^-^- p v , (163) 
where C = (2wa u ) 3j/ ' 2 e~ lEl ' T " . The integration over t leads to the expression 



T = Ce^u{p,) llx {l - 75 )w(p,)(0|J^_ A (0)|7r)e-^ 



sin(o;T/2) 



(164) 



The fraction in the bracket coincides with the delta function 2ir5(oj) of rep- 
resenting the energy conservation at the limit T — > oo and its deviation from 
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the delta function at a finite T is negligible, if its behavior near the roots 
of of oj — is normal of having a finite derivative. In this situation, the 
finite-size correction becomes negligible. The finite-size corrections of the 
integral or the average of probability 

n 2 



|T| 2 = \C\ 2 \u(p,) lfl (l - l5 )u(p u )(0\J^ A m^)\ 2 e-^ 5r 



2 sin(u;T/2) 



(165) 



are determined by the roots of u = and the behaviors of the u around the 
roots. 

The normal root satisfies 

5p = 0, 5E = 0, (166) 

and agrees to the solution of the plane waves, u is expressed by the momen- 
tum of the pion and muon as 

cu = E w - £ M (p M ) - E v (p w - p M ), (167) 

The roots are on the ellipse of q '— p n — p^ 

x 2 y 2 if ( in \ 2 1 

- + - = i, I = 9l --|i-^j « = <,„, 

where the direction of p n is chosen to the x-axis. The derivative of the angular 
velocity with respect to the momentum is expressed with the velocities 

which has the magnitude 

(d \ 2 
g-ouj =(v u -v fl cos9) 2 + vl(l-cos 2 9), (170) 

where v^, and v v are velocities of the muon and neutrino and 6 is the angle 
between p n and p^. The u varies from u = to ou(p + dsn norma i) in the 
normal direction and the value is computed as 



<q 2 1 11 Id 1 
u(q + dsn normal ) = ds\/-^ — (l - -) + - > ds-, ^-u > -. (171) 

71 



Since the slope is not small, the integral of \T\ along the normal direction 
with the parameter s converges fast at the large T. Hence the ordinary pre- 
scription for the square of the fraction 

2^')'=MH (172) 

is valid. The finite-size correction vanishes. 

Next we study the roots of u = of Sp ^ 0. The u is expressed in the 
following form, 

2 

777, 

u = E n - E(p„) - v v ■ (p, - p») - -f . (173) 
The condition u = is fulfilled by the momentum q = p n — p^ on the ellipse 

(q x + v v (E n - k)) 2 (1 - vl) + (q y - (p^yf = ^ V " El [ - ^lA _ _ m ^ 

(174) 

where the x-axis is chosen in the direction of p v . This ellipse has a large 
major axis that is inversely proportional to (^) 2 , and a minor axis that is 
proportional to (— ). 

On this curve, there exists a special solution of, bp = 0, SE = 0, which 
has the finite derivative of oo along the normal direction as was computed in 
the previous part. There exist other roots of Sp ^ 0, SE ^ 0, away from the 
above root, which have different properties. Because the ellipse is extremely 
long, the derivative of u along the normal direction in the second solution 
becomes 

%-% <™> 

This value is extremely small and |T| 2 varies extremely slowly in the normal 
direction and the integral over s deviates from that obtained from the infinite 

2 

region by an amount that is determined by the ^T. This gives the finite-size 
correction in the region T < — The finite-size correction has the origin in 
the root of u = at SE ^ 0, Sp ^ 0, hence does not hold properties that are 
derived from the energy and momentum conservation. Thus the finite-size 
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correction is computed uniquely with S[T}. Having tiny mass, the finite-size 
correction for the neutrino appears in the macroscopic area. 

For the detection of the muon, on the other hand, the muon velocity 
replaces the neutrino velocity v v in Eq. fjl64j) and the ellipse for the muon 
obtained from u' — with 

u' = E n -E fl {p fl )-E u {p u )-v fl {5p), (176) 

does not have large major axis. Hence u' has a finite slope on the normal 
direction to the ellipse and varies steeply at the large T. The integral of \T\ 2 
over s converges fast and the formula of Eq. ( 1172)) can be applied. 

Thus the reason why there is the large finite-size correction in the proba- 
bility of detecting the neutrino was elucidated. Nevertheless, it is not easy to 
integrate the muon momentum and to obtain the finite-size correction using 
the expression Eq. (I163p . Especially Lorentz invariance is not clear in this 
method. Instead, the method of integrating the coordinates in |T| 2 at the 
end is suitable for a rigorous computation of the finite-size correction, which 
was applied in the text. 
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